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THE RELATIONSHIP- BETWEEN MICROSCALES 
AND WIND-WAVE SPECTRAL DEVELOPMENT 


Craig Emery Dorman 


Submitted to the Joint Committee in the Earth Sciences, 
Massachusetts Institute of Technology and Woods Hole 
Oceanographic Institution, in October, 1971 in partial 
fulfillment of the requirements for the 
degree of Doctor of Philosophy 


ABSTRACT 


The objective of this study was to describe the mechanics of wind wave 
generation and spectral development* Intermittency, high frequency micro¬ 
structure in wind and wave fields, and strong nonlinear coupling involving 
a wide range of scales are shown to be crucial elements in the transfer of 
momentum to, from, and within the wave field. None of these elements are 
included in available theories. 

Measurements of wave height and of the turbulent atmospheric and sub¬ 
surface boundary layers were made from a small surface following platform 
and from a stable 38.5m spar buoy. The structure of moving gust patterns 
(cat’s paws) is described and related to the generation of surface waves. 
Results from this and other background studies are then applied to a dis¬ 
cussion of spectral growth during a two day period of active wave 
generation. 

Cat’s paws contain ’bursts’ of intense turbulent stress and buoyancy 
fluctuations separated by quiescent ’intervals’. There is a difference of 
over three orders of magnitude in fluctuation strength between these 
features. Rapid growth rate generation of high frequency surface waves 
and atmospheric turbulence occurs during the bursts. The resultant micro¬ 
scale components aid the growth of lower frequency instabilities by strong 
nonlinear coupling between scales of motion and by acting as drag or 
roughness elements. Evidence of strong coupling between frequency bands 
and of weakly resonant capillary-gravity wave interactions is presented. 
Thermal stratification has a strong influence on fluctuation magnitude 
and can delay the onset of surface wave generation. 

Major spectral growth is highly unsteady. Much of the momentum flux 
from air to sea occurs during intermittent events that are similar in 
nature to cat’s paws, and goes directly into high frequency waves. The 
bursts occur predominantly over large groups of surface waves and involve 
strong nonlinear interactions between media and frequency bands. The 
long-term equilibrium balance between wind and water is disrupted by 
variations in surface currents. There are ’critical’ wind speeds charac¬ 
terized by anomalous relationships between parameters of the predominantly 
logarithmic velocity profile. 



f 





r: 


«4ri 



5. 


I. INTRODUCTION 


Description of the Study 

This paper describes an experimental study of the mechanics of wind- 
wave generation and the spectral development of sea surface waves. Em¬ 
phasis is placed on the time history of generation and development, and on 
the nature and time scales of air-sea coupling. The basic objective was 
to describe and illustrate the physics involved in the transfer of 
momentum to, from, and within the wave field. The underlying hypothesis 
is that this process is highly nonlinear and depends strongly on the 
microstructure of the wave and wind fields. The experiments were con¬ 
sequently focused upon investigation of intermittent interactions between 
variables and scales of motion. 

Two separate experiments are described. The first is a detailed 
investigation of the structure of moving gusts. The second is an anal¬ 
ysis of spectral growth during a two day period of increasing wind. 

Both separately and together their results indicate that wind-wave gen¬ 
eration and spectral development involve strong nonlinear interactions 
with the wind and current fields; representation of these processes in 
terms of linear or weakly nonlinear coupling is inadequate. 

Background 

The available linear and weakly nonlinear wave generation models 
fail to yield satisfactory estimates of wave growth. In addition there 
exist unexplained anomalies such as spectral overshoot and the degener¬ 
ative effect of surface slicks on long waves. This combination of qual- 
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itative and quantitative unpredictability indicates that the physics of 
wave generation and spectral development is still inadequately under¬ 
stood. Until the correct physics is described, there can be little sig¬ 
nificant progress in either theoretical modeling or forecasting. 

, Since 1957, most wind-wave generation experiments have been designed 
around the linear Miles-Phillips models (Miles [1957], Phillips [1957, 
1966]). Early studies compared observed spectral growth with theoretical 
values; they soon showed that spectral shape and initial growth rates 
were reasonably well predicted, but that growth rates of the major por¬ 
tion of the spectrum were severely underestimated. Consequently, more 
recent work in the low frequency range has concentrated on the direct 
relationship between wind and wave fields in attempts to determine the 
adequacy of the theory’s representation of momentum transfer mechan¬ 
isms. The high frequency range (f >Ihz) has been treated as a separate 
entity, primarily relevant to initial generation. 

Work based on the simplified linear model has essentially dealt 
with the general or time-averaged characteristics of wind-wave coupling 
and spectral development. It has neglected the relationship between 
portions of the wave spectrum, and in particular has deemphasized the role 
of high frequency wavelets and the coiplLng between waves and turbulence. 
There has thus been very little work encompassing the full range of spec¬ 
tral growth and relating it to interactions within and among the wind, 
wave, and surface current fields. Hasselmann’s Hamiltonian type model 
(1967) was so general (and unfamiliar to field workers) that it also 
failed to inspire experimental investigation along these lines (with 
the exception of some ongoing parts of JONSWAP). 

A second experimental approach to the problem of spectral growth 
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has been through parameterization of wind and wave fields and' the devel¬ 
opment of spectral transfer functions for use in forecasting. Such studies 
have described the generally logarithmic nature of the atmospheric 
boundary layer over water, specified the characteristics of wave spectra 
in some detail, indicated the significance of energy transfer between 
wave components, hinted at the importance of high frequency elements, 
and yielded partially valid prediction schemes. They provide little 
direct information about the physics of momentum transfer processes, but 
rely upon theory and experience for guidance in the choice of empirical¬ 
ly applied parameters. 

It is highly unlikely that either forseeable modifications to the 
available models or additional parametric type studies will correct the 
known theoretical deficiences. Further advances will require a more 
complete understanding of the physics of generation. Important elements 
in a reevaluation of the wave ^noration process must include the detailed 
structure of the turbulent boundary layers and nonlinear couplings, both 
weak and strong, within and among the three fields of motion. Theoreti¬ 
cal and experimental studies have shown that such couplings do exist. 
Turbulence is known to be intermittent and strongly dependent on boun¬ 
dary conditions; a large number of scales are involved in the dynamics of 
flow at typical geophysical Reynolds numbers; and the generation of waves 
is inseparable from the generation of turbulence. These considerations 
imply that wave and turbulence generation may be a strongly nonlinear pro¬ 
cess involving a large variety of forces and scales, and that the evalu¬ 
ation of experimental results in terms of mean values and long averaging 
times can easily obscure the physics. Viewed in these terms the time his¬ 
tory of events during wave generation becomes an important factor; the 
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microstructure of the air and water boundary layers and the wave field 
reassumes great significance; and the ability to distinguish the differ¬ 
ences between generation, transfer, and detritus becomes a major experi¬ 
mental goal. 

The pair of experiments described herein was designed to investigate 
the waves and near-naviface boundary layers and to determine the nature 
and time scales of the important processes in wave generation and spectral 
development. Events ranging in duration from less than one second to 
two days were examined to obtain answers to the following questions: 

1. What are the important time scales in the transfer of 
momentum to, from and among surface waves? 

2. What are the nature and magnitude of the fluctuations 
involved in these transfers? 

3. What is the nature of the coupling between media, and 
between scales of motion both within and between media? 

4. What is the * imprint’ of wave generation in the air and 
in the water; and can this imprint be used to distinguish 
events of importance in the generation process (i.e., both 
to predict and to permit conditional sampling)? 

5. Is there a similarity between scales, and if so, for what 
variables and over what range of scales does the similarity 
exist? 

The first experiment was designed to isolate and intensively analyze 
an event of known significance - the ’cat’s paw’. The cat’s paw is an 
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isolated parcel of atmospheric turbulence. It contains anomalously large 
velocity and temperature fluctuations, and can be visually delineated by 
the *dark’, moving patch of high frequency waves which it creates. The 
second experiment yielded data from a larger system. It was designed to 
extend the time scales and number of forces involved to values typical 
of open ocean spectral development. The experiments are first individually 
discussed in separate papers; then the above questions are answered and 
the development of a wave field is described in terms of the results. An 
outline of experimental details is contained in appendices; complete des¬ 
criptions and all significant data are available in additional reports of 
the M.I.T. Buoy Project. 
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II, THE STRUCTURE OF MOVING GUST PATTERNS (CAT’S PAWS) 
AND THEIR ROLE IN WAVE GENERATION 


Abstract 


Measurements of wave height, turbulent velocity, and temperature 
fluctuations within 1 meter of the sea surface were taken from a small 
surface following platform. The resultant data show that moving gusts 
are highly structured, consisting of bursts of turbulent stress separ¬ 
ated by relatively quiescent intervals. Magnitudes of stress and*buoy¬ 
ancy fluctuations vary between bursts and intervals by over 3 orders of 
magnitude. The bursts are associated with the generation of very high 
frequency wavelets which in turn serve to couple the wave field ^to the 
mean air flow, A large range of scales of motion interact, leading to 
a lognormal distribution of turbulent velocity fluctuations and to ener¬ 
gy transfer among wave components. 


Introduction 

Gust patterns moving over the sea surface (cat’s paws) are well- 
known to sailors, both those seeking a strong wind and those trying not 
to capsize. Cat’s paws can be located by their imprint on the water 
surface, namely the change in surface reflectivity caused by the high fre¬ 
quency wavelets generated. Associated with these gusts are large magni¬ 
tude, high frequency, turbulent stress and buoyancy fluctuations. Cat’s 
paws are patches of intense air-sea interaction with correspondingly large 
transfers of momentum, heat, and moisture. This study reports an attempt 
to measure, describe, and find relationships among the fluctuations with¬ 
in a cat’s paw, with particular reference to the mechanisms involved in 
momentum transfer and wave growth. 

Incorporated in a cat’s paw are three elements of wave generation 
and spectral development that are neglected in the available modern models 
(Jeffreys [1925, 26], Phillips [1957, 1966], Miles [.1957], Hasselmann 
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[1967]). These are intermittency of generation,' the role of wavelets in 
spectral development, and direct, strong nonlinear coupling between scales 
of motion and between the sea surface and atmospheric turbulence. The cat’s 
paw is itself an intermittent or burst type phenomenon; the activity with¬ 
in it is qualitatively and quantitatively different than that observed 
on the average. As an experimental objective, a gust serves as a realiza¬ 
tion one type of wave generation event, with constraints determined by the 
environment, at Reynolds numbers unattainable in the laboratory. 

Although the above three elements have been largely neglected, sub¬ 
stantial eveidence of their significance is available in the literature. 

Munk noted in 1955 that while high frequency waves account for an insig¬ 
nificant portion of total wave energy, they play a disporportionate role 
in both wind stress and the development of the spectrum. This had been 
exemplified by Van Dorn’s (1953) experiment, where the presence of a sur¬ 
face slick significantly reduced setup in an artificial pond. This 
effect of slicks was recently reemphasized (Ruggles [1969], Barger, et al, 
[1970]) by a demonstration that removal of wavelets resulted in a re¬ 
duction of wave energy accompanied by an increase in mean wind, decrease 
of roughness length, and significant reduction of coherence between 
wind and waves. 

Stewart (1961, 1967) has long maintained that much if not most of 
the momentum transfer from air to water is effected by wave drag. Van 
Dom (1953) and Wu (1969 a, 1969 b, 1970) attribute this wave or form 
drag to the high frequency wavelets superimposed on larger waves. Wu 
further relates the wavelets to roughness height and states that the mo¬ 
mentum flux to wavelets accounts for a substantial fraction of the wind 
stress; the stress appears as a normal component to the wavelets, but 
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feels like a tangential stress to the underlying waves. Longuet-Higgins 
(1969a) has shown that such a fluctuating tangential stress, occurring 
selectively near the crests and leading faces of longer waves where the 
shorter ones are ^trapped^ (Munk [1955], Cox [1958]), is dynamically equi¬ 
valent to a normal pressure fluctuation in quadrature with them. 

The role of turbulence and turbulent Reynolds stress in wave gen¬ 
eration is commonly viewed as minor or negligible. Miles-Phillips theory 
holds that the waves do not affect turbulence; Hasselmann includes a weak 
interaction through a disturbance of mean flow by one wave component, 
interaction of the disturbance with a turbulence component, and a feed¬ 
back of the resulting pressure perturbation into another wave element. 
Stewart (1967) has shown that wave energy is typically comparable with 
the energy of turbulence in the entire troposphere, and consequently has 
suggested that whxle turbulence can serve as a catalyst for transfer of 
momentum from the mean flow, it can by itself add very little to the wave 
energy. 

Two of the more important properties of atmospheric turbulence with 
respect to wave generation are intermittency and the associated large 
growth rates of localized high frequency instabilities (Greenspan and 
Benney [1963]). Laboratory experiments (e.g., Klebanoff, et al; [1962], 
Kline, et al, [1967]) have shown that the progression of events leading 
to instability and transition of laminar boundary layers involves suc¬ 
cessively the formation of travelling Tollmien-SchlichtSng waves, develop¬ 
ment of spanwise disturbances leading to formation of local shear layers, 
and explosive development of the flow into short turbulent length scales 
in bursts of high frequency instability fluctuations. The turbulent 
boundary layer is itself intermittent, with periods of intense activity 
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separated by relatively quiescent intervals (e.g., Frenkiel & Klebanoff 
[1971], Rao and Narayanan [1971]). Reverse transition is similarly 
abrupt (Steiner [1971]) and is due to the restoring action of buoyancy 
forces. Furthermore, Klebanoff, et al, (1962) demonstrated that the 
spatial structure of turbulent breakdown is stroiigjy dependent on boundary 
irregularities, and Guptai et al, (1968) showed that unstable oscilla¬ 
tions irdiced by the simultaneous introduction of Tollmien-Schlichting 
waves and free surface waves resulted in a large, unpredicted amplification 
rate of small wavelengths. 

The combined effects of intermittency, high growth rates, and boun¬ 
dary influence were illustrated (Mollo-Christensen [1970]) by filtering, 
squaring, and plotting selected frequency components from wave gauge and 
anemometer signals recorded during a period of wave generation. The plots 
showed that ripples and high turbulent wind energy were associated with 
groups of large amplitude, low frequency waves. This suggests that 
generation of waves and turbulence are coupled, intermittent, and asso¬ 
ciated with the presence of wavelets on longer wave crests. The basic 
roles of the high frequency waves are the coupling of wind and wave 
fields and the enhancement of low frequency instabilities; the wavelets 
are generated by bursts of high frequency turbulence, and growth rates 
in these local instabilities are extremely rapid and due to strong non¬ 
linear coupling. In turn the turbulence patterns are controlled by the 
major wave field. The combined result is a strongly coupled system with 
growth rates much larger than those predicted by linear or weak coupling 
theory. 

Implicit in this outline of spectral development is the participa¬ 
tion of a large range of scales, coupled by nonlinear interactions. The 
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high frequency waves are in the capillary-gravity and capillary wave 
ranges. As pointed out by Stewart and Manton (1971), the important 
time scale for generation by turbulent pressure fluctuations (Phillips* 
resonance mechanism) is that of the wave group generated, i.e., the wave 
frequency. Similarly, the large amplification rates in shear layer in¬ 
stability (Greenspan & Benney,[1963]) have a typical time scale of a 
half period of the primary wave. Water surface waves of the capillary- 
gravity type have second order resonant interactions (McGoldrick [1965], 
Simmons [1969]) with effective time scales of two to three waves periods. 
Turbulent bursts or spots are similarly of short duration. It can 
therefore be anticipated that the generation of wavelets is a short¬ 
lived but intensive phenomenon which is capable of rapidly depleting the 
locally available energy of the mean flow. At the same time, as pointed 
out by Mollo-Chriotensen (1971), **work done on the water to make waves must 
equal a stress times a strain, and the reactions of the stresses on the 
water act on the air to produce turbulence.** 

The intermittency and selectivity of the location of high frequency 
waves on the underlying long ones is only partially attributable to the 
coupling with atmospheric turbulence.' Variations in surface tension near 
the crests of steep waves produce a local increase in normal stress which 
in turn generates capillary waves (Longuet-Higgins,[1963], Crapper [1970]). 
The short waves steepen on the forward faces of longer waves due to hor¬ 
izontal contraction; this leads to breaking and the transfer of mass and 
momentum from the short to the longer waves (Longuet-Higgins [1969b]). 
Hasselmann (1971) has shown that the work done on the long waves by the 
**interaction stresses** is balanced to first order by the loss of poten¬ 
tial energy involved in mass transfer, and to second order results in 
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a slight attenuation of the longer waves. The basic physical idea sug¬ 
gested by Longuet-Higgins has not been negated by Hasselmann’s rescaling 
of the terms; there does exist a considerable energy interchange between 
scales of motion within the water, and this interchange is strongly depen¬ 
dent on the presence and location of high frequency waves. 

A similar interplay occurs in the atmosphere. Mollo-Christensen 
(1970, 1971) has shown that when the velocity field is described in terms 
of three components - mean flow, long wavelength and short wavelength 
fluctuations - the resultant energy equations involve interaction or pro¬ 
duction terms, in the form of stress times strain, among all three com¬ 
ponents. If the interactions occur with the proper phase relationship, 
a primary instability can be enhanced by the work done by Reynolds' 
stresses of the secondary instability, and both instabilities can for 
a while gain energy from the mean flow. This type of interaction can 
occur similataneously on many scales, in both the wind and the wave 
fields. The essential points are that Reynolds stresses are inseparable 
from the generation of fluctuations, and that by nonlinear couplings between 
scales, the growth of short waves .can enhance the growth of longer ones. 

On the basis of these considerations from the literature, it is 
suggested that the major reason for the inability of available theories 
to satisfactorily predict spectral growth is that they fail to include 
the correct physics of the momentum transfer process. The inclusion of 
strong nonlinear couplings tremendously increases the complexity of the 
problem ( and makes it as yet mathematically intractable), but appears 
to be necessary for a proper understanding of wave generation. 

The cat's paw is a readily accessible phenomenon that provides a 
realistic geophysical framework within which to initially investigate these 
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ideas. The wind tunnel is too small for studies at the Reynolds num¬ 
bers and hence the range of scales involved in the proposed mechanisms, 
and the dynamic range of wave height on the open ocean makes detailed 
near-naviface measurements exceedingly difficult. This experiment was 
designed to investigate the time history of cat’s paws and to determine 
the scales involved in them. The primary features sought were the rel¬ 
ative magnitudes of the fluctuations in air and water and an indication 
of coupling between media. Thus the variables of major concern were 
turbulent Reynolds stress and buoyancy fluctuations, and high frequency 
wave activity. Interactions between scales of motion were examined in 
terms of the basic cat’s paw structure, turbulent velocity fluctuations, 
and nonlinear wave interactions. Specifically, answers were sought to the 
following questions: 

Is there a demonstrable relationship between the generation 
of turbulence and the generation of high frequency waves? Is the 
process intermittent or is it typically near the mean or average 
flux? 


What is the evidence for coupling between scales of motion? 


On what time scales does instability-associated high growth 
rate generation occur? 


What is the imprint of a generation ’event’ in the air and 


water? 
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The Experiment 

The measurements were taken from a small (1 m long) surface following 
catamaran (cat; Figure 1). The cat was free in pitch, roll, and heave 
(design natural frequency 1.5 hz) and restrained in yaw and horizontal 
position. It could be oriented into the wind through the full range of 
acceptable fetch directions. The experiment was conducted in March 1971 
from the dock of the M.I.T. Yacht Club on the Charles River, Cambridge, 
Massachusetts. 

Instruments included 6, DISA type 55A81 hot film anemometers, a capac¬ 
itance wave gauge with .076 cm diameter sensor, and VECO 42A70 (air) 
and Fenwall K496 (water) thermistors balanced against precision resis¬ 
tors in a differential bridge circuit. Subsidiary equipment included 
an accelerometer for sensing catamaran motion and cup and propeller anem- 
ometeis for determining mean wind. Data were recorded on magnetic tape 
in analog mode, accompanied by a voice narration that was later used 
for subjective conditional digitization. 

Two of the hot films were orthogonally mounted to permit determina¬ 
tion of horizontal ( U ) and vertical (W ) velocities and hence Reynolds 
stress (-/^o'W ). The sensors were attached to the forward vertical cat 
support to minimize tilt (<^1® mean, <“3° fluctuating). Table 1 lists the 
instrumental arrangement and mean conditions for the seven 2-hour runs 
that were analyzed. Data from runs 2, 4, and 6b, which represent the 
full range of conditions encountered, were selected for presentation 
herein. 

The data were analyzed using a combination of digital and analog 
methods. Pre-analysis, demodulation, and filtering (100 hz low pass) 
were analog. The data were then digitized at 250 hz and recorded on 7 
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Figure 1 


Catamaran and Support Structure 
























TABLE 1: Catamaran Runs (2 hour duration) 
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track digital magnetic tape. Final analysis was accomplished using lo¬ 
cally generated software on the M.I.T. IBM 360. All available channels 
were analyzed simultaneously for each run segment of interest. Exten¬ 
sive use was made of graphical output from a Stromberg-Carlson 4020 
plotter. 

The major contraints on the analysis were the cat motion and the 
nature of the cat’s paw itself. Cat motion could be successfully removed 
from the data, but the very fact that it did follow the low frequency 
waves prevented meaningful quantitative estimates of spectral growth. 

The cat’s paw will be shown to consist of patches of strong fluctuations 
separated by relatively quiescent intervals. Thus time averaging and 
the selection of data segments for analysis had to be very carefully 
done to avoid distortion of fluctuation magnitudes and to prevent ob¬ 
scuring events of importance. The time history of events proved to be 
of crucial importance. Further, the momentum transfer processes involve 
both internal and external interactions; they are highly non-stationary, 
non-homogeneous, and anisotropic, so that the standard statistical theo¬ 
ries of turbulence are inapplicable and time averaging could not be sub¬ 
stituted for ensemble averaging. Therefore data from any one gust pos¬ 
sesses little if any statistical significance, and conclusions could be 
drawn only on the basis of the change with time in the individual prop¬ 
erties of a variable and in the relationship between variables, and then 
tested on the basis of similarity among a number of identically treated 
events. Each cat’s paw, and each segment within a cat’s paw, must be 
viewed as a single realization. No meaningful error bars or confidence 
limits can yet be placed on the data since a statistically insufficient 
number of events has been evaluated. The data were subjectively viewed 
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in terms of the relative changes with time, the sequence or pattern of 
events, and the relationship between variables in both time and frequency 
domains. Reported magnitudes can be considered only as order of magnitude 
estimates. This is particularly true of wave*autospectra, where filtering 
and the necessarily very short averaging times preclude even rough estimates 
of the magnitude of spectral growth and decay. 

Emphasis was placed on understanding the intermittent nature of tur¬ 
bulent fluctuations, their interaction with the free surface, and changes 
in the presence and structure of capillary-gravity and pure capillary 
waves. These features were then evaluated in terms of their significance 
to the generation of surface waves. 


Results and Discussion 


A. Cat*s Paw Structure in Time and Frequency . Separating velocity 
and temperature fields into steady (U,T) and fluctuating (u’,T’) components, 
the equation for the total kinetic energy of turbulence can be written 
(Lumley and Panofsky [1964]) 


® d) © 

where Einstein’s summation convention holds, yU is dynamic viscosity, and 
To the temperature of the undisturbed atmosphere. is commonly called 

the production term; it represents the turbulent energy produced by fluc¬ 
tuations detracting energy from the mean motion by instability mechanisms 
as work=stress x strain. represents work done against buoyancy in 

terms of heat flux, ^p '^'^i ’ upward heat flux acts as a source, 

a downward flux a sink of kinetic energy. © is the rate of working 
against viscous forces and represents dissipation of energy to heat (£ ). 
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The term from A is the Reynolds stress ( r ); it commonly 

represents momentum flux from the air to the sea surface (downward flux 
for > 0). This transfer across the boundary acts as an additional sink 
term as seen from the viewpoint of atmospheric turbulence. Similarly, 
it is apparent that the reaction stresses caused by the generation of 
waves acts as a source term for the generation of turbulence. Thus gen¬ 
eration of waves and generation of turbulence are inseparable. 

Horizontal (U) and vertical (W) velocities, termperature (T), wave 
height ), Reynolds stress (z“ = -/>u’w’), and buoyancy • (/^ = w*T’) were 

measured or computed and related in time and frequency domains to depict 
cat’s paw structure. All data were taken at a single horizontal location, 
and w’ and T were measured at only one height in air and water; it was 
therefore impossible to determine a meaningful quantitative energy bal¬ 
ance, although the qualitative interplay of production and stability terms 
was apparent. 


A.l. Individual Variables . Figures 2 through 7 (Figs. 2 - 
7) depict the imprint of three gusts. Fig. 2, from run 2, shows the 
fluctuations In ^ ,11,1^7:, and /S for 40 seconds (4, 10-second time per¬ 
iods [T.P.]) starting just in advance of the leading edge of a gust and 
ending aft of the trailing adge. Fig. 3 shows samples of the computer 
generated (SC-4020) plots from which Fig. 2 was drawn (W is also shown). 

(The arrows and slant numerals on the t plot of Fig. 2 indicate the 2- 
second and 1-second segments chosen for detailed analysis.) Fig. 4 and 

Fig. 5 are similar representations of a cat’s paw from run 4. ^ >U, W, X , 

2 

a switch or intermittency function of -c 100 dynes/cm , T, 
r ^ are shown. The three non-consecutive T.P. are from before the 
gust, near the leading edge, and about 5 seconds later, still within the 
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gust. Fig. 6 and Fig. 7 are similar diagrams from run 6B. In run 6B 
there are three consecutive 10-second T.P. starting just at the leading 
edge, and one T.P. from 10 seconds after the gust. These gusts were 
selected to illustrate the full range of velocity and bulk stability 
values encountered (Table 1). 

The basic features in all three runs areas follows. The leading 
edge of a gust is marked by rapid increases in both mean and turbulent 
wind velocities. Within a very short time (< 1 second), the air temper¬ 
ature drops as the boundary layer is mixed and destabilized by the tur¬ 
bulent stresses. T continues to drop as the winds increase; it starts 
to increase after U has reached its maximum value. There is commonly 
a secondary increase in U near the trailing edge that causes an additional 
drop of T. In run 6B, this second T drop is more rapid than the first, 
since the boundary layer has not yet recovered its initial static sta¬ 
bility and the stress fluctuations are consequently more effective in 
mixing. After gust passage, T rapidly rises to its initial value as the 
boundary layer restabilizes. There is no clear-cut relationship be¬ 
tween atmospheric variables and the low frequency wave field or water 
temperature. 

Both *2" andyS traces are highly intermittent. There are bursts of 
high frequency, large magnitude fluctuations, separated by quieter 
intervals. Further, there are well-defined sub-bursts and intervals 
within the larger bursts. The r fluctuations are a maximum during in¬ 
creasing wind, but are not restricted to these periods. They are strong¬ 
est and contain most high frequency components not at the gust leading 
edge, but rather in the central bursts, after the boundary layer has 
been mixed and near the point of U . Bursts in r are generally asso- 
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ciated with similar but shorter and weaker events in /e ; such 
coupled events are followed by T decreases. Bursts inalso occur 
independently or with relatively stronger magnitude than the fluctua¬ 
tions. These events are followed by a levelling off or rise in T. The 
stress dominated fluctuations destabilize and mix the boundary layer, and 
buoyancy forces act rapidly to restratify it. The buoyancy forces are 
relatively stronger during run 6B, with high bulk stability. Here the 
sub-bursts are clearly evident. Comparison of the three runs, and of 
the two T.P. from run 6B, indicates that there are both stronger and 
higher frequency ^ andfluctuations during higher winds. In the stable 
run 6B, ^ fluctuations die out before those in ^ ; in nearly neutral 
run 2, even small z: anomalies cause variations in ^ water fluctu¬ 

ations are also seen to be intermittent. They appear to occur during 

C\ 

atmospheric bursts and under breaking waves, "^^^ter 

very small, and strongly coupled to the wave field (i.e., they are due 
to motion of the thermistor through the vertical stratification); a 
small mean decrease results from mixing during the passage of 

the gust. (The spikes on the traces of Fig. 5 are due to 

least bit noise in the digitizer.) 

The overall pattern indicates that fluctuations in the leading edge 
of the gust cause initial destabilization of the stratified boundary 
layer. This is followed by intermittent large magnitude instability 
fluctuations.‘ Buoyancy forces start to restabilize the boundary layer 
during quiescent intervals, and fully restore stability within 10-15 
seconds after the gust has passed. 

2 

The stress mean value, are on the order-of 10^ dynes/cm 
for run 2, 5 x 10® for run 4, and 10® for run 6B. is approximately 
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^ 1 Q ^ cm- ®C/sec for run 2, and +10 ^ for runs 4 and 6 B. These values 

are consistent with the differences in mean conditions and comparable 

with the commonly assumed values used in dynamic modeling of the exchange 

processes (r = 1 dyne/cm^, H = 3x10 gm-cal/cm-sec). The ’instantaneous’ 

3 2 

(250 hz) values commonly approach or exceed 10 dynes/cm and 

2 

3x10 cm-®C/sec in the bursts. The strong high frequency fluctuations 

thus have amplitudes typically three orders of magnitude greater than 

the mean, or departures from the mean of 10 to 50 standard deviations. 

Similar types of measurements taken 5.6 to 22.6 m over a wheat field 

during unstable stratification (20 hz sampling rate) (Haugen, et al, [1971]) 

2 2 2 

yielded burst magnitudes of r - 10 dynes/cm ,/^ -- 3x10 cm ®C/sec. Burst 
duration was roughly 10 seconds, and < U >^ 7 ;; 600 cm/sec. This data also 
showed the presence of bursts within bursts. McBean (1970) also found 
similar values over land, with a strong effect of stability (momentum 
transfer efficiency decreased and heat transfer efficiency increased with 
decreasing stability). The major qualitative patterns of intermittent 
gusts are thus similar over land and water. Quantitatively, the time 
interval selected for averaging is deterministic in estimates of instan¬ 
taneous fluctuation magnitude (Section C). A 15-minute averaging period 
for the cat’s paw data would have increased the above estimates by at 
least an order of magnitude. The difference in magnitude between cat 
and land results is an indication of the significaice of stability and 
the enhancement of frequency and magnitude of stress fluctuations caused 
by the participation of a fluctuating boundary capable of supporting 
well-defined high frequency waves. 

Fig. 8 shows some typical histograms of the individual variables 
from the 10 second T.P. of all three runs. Fluctuation magnitudes and 
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SAMPLE HISTOGRAMS FROM RUNS 2,4, a 6B 

10 SECOND SEGMENTS 

Distributions in 7o of total samples (2520) per interval 

_ Figure 8. _ 
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variance are much larger during tha,n before or after the gusts; the peak¬ 
edness of the ^ y /8 , and w histograms is due to the high frequency 

fluctuations and consequent large number of axis crossings. U and T’ 
distributions are more unimodal after than during the gust. There are 
few other basic differences between the T.P.'s for any one variable when 
the relatively long lO-second intervals are used. 

Fig. 9 shows typical autospectra based on the 10 second T.P. (6 degrees 
of freedom; 90% confidence limits are .63 x 10^ ^ to 1.74 x 10^ for a spec¬ 
tral density estimate of 10^). With this comparatively long averaging 
time there are only generalized differences in frequency content between 
gusts and mean conditions. With the exception of the wave gauge data, 
the spectra are fairly white, with (note run 6B data) a slightly more 
rapid drop-off with increasing frequency in the T.P. after the gust. As 
also shown in the time plots, the high frequencies are augmented during 
the bursts. The differences in magnitude between gust and no gust con¬ 
ditions are apparent; the peaks of (f) and (f) are 3 decades 

lower after the gust. This basic similarity of pattern, but difference 
in high frequency content and major difference in intensity reflects the 
patterns displayed in Fig. 7 where the T and /0 scales have been deter¬ 
mined by the computer on the basis of the 10-second extrema. Wave propa¬ 
gation and the presence of long period wave groups prohibit any deductions 
from magnitudes of the wave spectra. 

Fig. 10 displays the changes in the 0.5 second vertical profiles of 
U from run 6B. It also shows U ^ 0,5 sec and a measure of the turbu¬ 
lent kinetic energy production term. Profile perturbations, the devel¬ 
opment of unstable shear layers, and turbulent production occur predom¬ 
inantly during bursts. As with the r and^ fluctuations, production has 
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Figure 10 
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a range of i 3 to 4 orders of magnitude and is highly intermittent. The 
comparatively long averaging times CO*5 sec) are incapable of showing 
the extremely rapid fluctuations within the bursts; instead, the profile 
changes are similar to those that would be produced by a set of eddies 
propagating past the sensor array. Again, the maximum turbulence pro¬ 
duction occurs after the leading edge in T.P. 2, i.e., after the atmos¬ 
pheric boundary layer has been destabilized. The profile remains nearly 
steady, and production nearly zero, after the gust has passed and thermal 
stratification has been reestablished. Both stability and intermittency 
are of major importance in gust time history. 

A sample time history of three cat’s paws thus indicates that gusts 
are highly structured. Averaging techniques such as spectra and histo¬ 
grams tend to obscure the significance of the microstructure. and yS 

fluctuations in bursts are at least 3 orders of magnitude above the 
mean. The localized instabilities also contain a larger range of fre¬ 
quency components than is the case in quieter periods and intervals. 
Production of turbulence and mixing in the boundary layer are strongly 
associated with the intermittent bursts; r fluctuations are largest 
after the mixing and destabilization associated with the leading edge. 
Participation of a large range of scales of motion and of the fluctuating 
water surface enhances the burst magnitudes. The use of long averaging 
periods obscures the details of the transfer processes that are dependent 
upon these strong localized instabilities. 

A.2. Relationships between Variables . Relationships were 
depicted by joint histograms (JH) and time plots in the time domain, and 
by coherence (Coh) plots in the frequency domain. This duality of re¬ 
presentation helped to lend confidence to conclusions based on a limited 
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amount of data. Coherence may be defined as 

-i. 


It is represented in terms of its magnitude and phase, 

(rt = Y,yO . 

It was computed using fast fourier transform techniques. In the conven¬ 
tions adopted, means leads . Averaging was over 25 adjacent 

frequency bands for 10 second T.P., 9 for the shorter intervals.* As shown 
by McBean (1970), the’spectral correlation coefficient’ (or Coh) is an 
actual measure of transfer efficiency as a function of frequency (or scale). 


a. 10 Second Intervals . Fig. 11 shows JH(U,T’) 

(with 10 second mean and trend both removed from T) from run 6B. Asso¬ 
ciated with the decrease and increase of T > 10 sec is the boundary layer 
destabilization and mixing by the leading edge, violent mixing during 
the major bursts, and subsequent restabilization of a nearly constant temp¬ 
erature after gust passage. The top of Fig. 12 shows a similar pattern 
in the frequency domain using Coh^.',^ from run 4. Before the gust the 
waves and temperature fluctuations (above 5 hz) are quite coherent and 
nearly out of phse. This is the pattern to be expected if the catamaran 
rides up and down through a thermally stratified layer that follows the 
pattern of the waves. (Ramzy and Young [1968]). Near the leading edge 

is slightly lower and extends to lower f requencies ; begins to 
shift, particularly in higher frequencies, as the atmosphere is mixed. 

Later in the gust the temperature field has been effectively inverted, and 
T and yj are more nearly in phase. The lower half of Fig. 12 shows Coh^^ 
from the same T.P. Here is much lower but shows effectively 
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the same temperature inversion. 

Fig. 13 and Fig. 14 show JKQt ^ ) and Coh^^ from run 6B; Fig. 15 

shows JK( r ,^ ) from run 2. These diagrams show that the mixing and 
tumbling is effected by' T fluctuations destabilizing the boundary 
layer, and that the following restabilization is due to the predominance 
of buoyancy forces. The change of the patterns with time is most clear 
in the more stable run 6B (note the change in scales); in the nearly 
neutral run 2 the fluctuations show less correlation, are more randomly 
scattered. Here the stress fluctuations predominate and remain intermit¬ 
tently large even after the gust has passed, while in run 6B the buoyancy 
forces rapidly debilitate the momentum flux as the gust departs. The 
coherence diagrams also indicate the phase reversals of stress and buoyan¬ 
cy fluctuations as the boundary layer is turned and mixed. The rapid, nearly 
linear shift of ©r/S with f in the 9.5-17 hz range of T.P. 4 implies that 
fluctuations of similar frequency are ’slipping’ past each other, or 
effectively that in the presence of strong stratification, weak z: fluc¬ 
tuations are unable to overturn and mix the stable boundary layer. 

Boundary layer stability is effective in retarding momentum flux 
by turbulent velocity fluctuations. This effect is particularly large 
near the leading edge of cat’s paws. Wave generation during high fre¬ 
quency instabilities will thus be less efficient during stable than 
neutral stratification, and delayed or reduced by the conversion from 
kinetic to potential energy involved in disruption of the stability 
structure. 

Fig. 16 and Fig. 17 show the JH and Coh represetations of the re¬ 
lationship between U at two heights during run 6B. As vzith the other 
variables, fluctuations are largest at and just after the leading edge. 
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with a collapse of the histogram after the gust has passed. Yyy is sig¬ 
nificant only at the DC level (T.P. 1-3) and the reversals of ®uoCD*C.) 
demonstrate the large variations in mean wind velocity shear that occur 
during the gust. The generally low level of and lack of JH patterns 

over even this short distance are indicative of strong vertical struc¬ 
ture and the anisotropic character of the high frequency intermittent 
fluctuations. 

Fig. 18 and Fig. 19 are JH( T ) from run 6B and Coh from 
during and after the gusts of runs 2 and 6B and from all three T.P. of 
run 4. It is apparent that there is no effective evidence of strong 
wind-wave coupling during the gust when viewed from this averaging period. 
The effective generation mechanisms are obscured by the ’detritus’ or 
by averaging over events of radically differnt character. If anything, 

“c and Y are more coherent in the 7-12 hz range (of runs 2 and 6B) 
after than during the gust; and in this frequency band the continuous 
linear drift is evidence that the fluctuations of similar frequen¬ 

cy slip by one another rather than strongly coupling. The turbulent 
stress-wave relationship that will be demonstrated has again been obscur¬ 
ed by considering averaging intervals much longer than the effective inter¬ 
mittent events. 

b. 2 and 1 Second Intervals . Two second segments 
containing primarily ’burst’ or primarily ’interval’ type fluctuations 
were selected from runs 2 and 6B. The location of these 2-second T.P. 
is shown in Fig. 2 and Fig. 6 ^ith the exception of T.P. 1 from run 6B 
which was taken about 3 seconds prior to arrival of the leading edge.) 

A portion of the results is shown in Fig. 20, JH( r , ^ ) from run 2, 
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and in Fig. 21, Coh^^ from run 6B. These figures show some overall 
similarity to those of the 10 second T.P.; some of the segments still 
contain both bursts and intervals, and this tends to partially blur the 
effects of the very high frequency fluctuations. 

Fig. 20 shows that the mixing and destabilization process is indeed 
associated with the bursts. Relatively weak, positively correlated fluc¬ 
tuations occur near the leading edge; mixing occurs during the next 
strong burst and is followed directly by negative correlation of r: 
and ^ and consequent further destabilization of the boundary. With the 
boundary layer inverted the buoyancy forces work toward restabilization 
and rapidly degrade the strength of the momentum fluxes. Both the mixing 
and the restabilization occur rapidly; the intial boundary layer struc¬ 
ture is wiped out by the leading edge but is never firmly reestablished 
during gust passage. 

The wind-wave coupling begins to be apparent on this time scale. 

Fig. 21 shows that stress fluctuations are intially (before the gust) out 
of phase with the waves (7-15 hz range). They are almost totally inco¬ 
herent during the leading edge, i.e., during destabilization, then 

increases slightly in T.P. 3. This T.P. contains several v/ell-defined sub¬ 
bursts and intervals so that the picture is confused, but is relative¬ 
ly constant at the peaks at 5.9 and 12 hz. During the major burst, 
the high frequency instabilities and high frequency waves serve to strong¬ 
ly and immediately couple the lower frequency fluctuations to the under¬ 
lying waves. The high frequency coupling is still not evident, but is 

high at and below the main wave period, with or the maximum 

stress occurring aft of the crests. In the immediately following inter¬ 
val the D.C. end drops off slightly and the phase shifts; in 
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the higher frequencies is aligKtly higher but the propogating cap¬ 
illary waves now ’slip’ with respect to stress fluctuations of the same 
frequency. 

In order to further decrease the blurring of burst and interval, 
segments approximately 1 second long were selected from within the 2- 
second T.P. (see Fig. 6). The results from run 6B are shown in Fig. 22- 
Fig. 24. T.P. 1 is prior to gust arrival; T.P. 2 just behind the leading 
edge; T.P. 3 during a short but intense fluctuation 2 seconds later; T.P. 

4 during the major burst and T.P. 5 from the quiescent interval 3-4 
seconds after T.P. 4. The pattern of events during this sequence may 
be described as follows. 

Prior to the arrival of the gust, t and ^ are nearly totally co¬ 
herent and in phase. X fluctuations are small, loosely coupled to the 
main wave field, and out of phase with the high frequency waves. With 
the arrival of the leading edge of the cat’s paw the stress fluctuations 
increase in frequency and magnitude and the boundary layer is tumbled; 
the T j jS fluctuations decrease in coherence and* become nearly out of 
phase; little effect is noted immediately on the water surface, as the 
turbulent kinetic energy is utilized to stir the boundary; waves and "r 
are incoherent. After the intial burst the air is mixed and high fre¬ 
quency waves are generated, associated with the short, strong bursts of 
momentum flux. T and ^ are again coherent, but out of phase, and the 
wind-wave coupling is apparently random. VTith the boundary layer de¬ 
stabilized, the magnitude of the stress fluctuations increases (T.P. 4) 
and the high frequency instabilities become strongly associated with the 
generation of surface waves. is quite low at frequencies above 15 hz 

with more variable than before but centered around O . 
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increases considerably in the ^ ^ZB hz range, but remains fairly 
low for lower frequencies. slips across part of the coherent ranges 

but shows a major reversal of trend in the center of the coherent band. 
During the quiescent period immediately after this event r fluctuations 
are much smaller and less strongly coupled to the short wave field, with 
slightly increased in lower frequencies, r j fluctuations are 

again strongly coupled, and out of phase due to the mixing of the preceed-- 
ing burst. This process is repeated during subsequent bursts, until 
the cat^s paw has passed the site and stability is reestablished. 

The change in character of the high frequency stress and wave pat¬ 
terns (Fig. 23) is strong evidence of the large growth rates associated 
with intermittent, large magnitude, high frequency instabilities. The 
generation of high frequency waves, as in T.P. 3 and especially in T.P. 4, 
is indeed strongly coupled to the generation of turbulence. The large 
value of ^ 17 ^ for lower frequencies in 2-second T.P. 4 indicates that 
the high frequency waves can rapidly and effectively couple the lower 
frequency waves and instabilities, and aid their-growth. It is obvious 
that little statistical confidence can be placed in low frequency measures 
based on such short time periods; neverthegss the change between time 
periods and averaging lengths, the similarity of this pattern in other 
bursts and gusts, and the enhanced high frequency 1-second T.P. 4 

lend strong subjective credence to the hypothesis. 

Figs. 23 and 24 also clearly show how wave generation is delayed by 
the stability structure. Generation processes act only after the boundary 
layer has been initially destabilized. The continuous Ozrj slipping in 
T.P. 5 (across the full width of the coherent band) appears to be charac¬ 
teristic of 'detritus*, in this case, of the propagation and advection of 
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the waves and turbulence generated during the high frequency, localized 
instabilities. The slipping pattern is similar in T.P. 4 except for the 
10 - 18 hz range, where levels off and reverses direction near the 

zero level. Most of the generation and momentum transfer occurs in this 
frequency range (see Fig. 27). 

Generation appears to be a very strong but selective and short¬ 
lived phenomenon. The time scales associated with the sub-bursts and 
instability associated transfers are very short, typically on the order 
of 1 to 1.5 seconds, which is also the typical period of the immediately 
underlying major wave field. The strongest coupling and high growth rates 
associated with the localized instabilities occur in a narrow high fre¬ 
quency band. Averaging times significantly longer than these scales can 
depict the general changes in atmospheric structure but fail to show 
the detailed nature of the coupling and transfer processes. The impor¬ 
tant scales are thus those of the waves, the bursts, and the internal 
fluctuations within the bursts. Momentum is transferred to the water 
at scales associated with these features. 

B. Interaction Between Scales . The microscales in both wind 

and wave fields account for only an insignificant portion of the total 
energy, yet as shown above they are the active participants in the local¬ 
ized high growth rate instabilities. They are effective in aiding the 
growth rate of lower frequency components and in coupling wind and wave 
fields, as was heuristically explained by Mollo-Christensen (1970, 1971) 
and as demonstrated in the time plots and coherence diagrams of the pre- 
ceeding section. Other evidence of nonlinear coupling between scales 
may be shown by an examination of the probability distribution of turbulent 
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velocity fluctuations and of the shape of the high frequency wave spectrum. 


B.l. Scale Interaction in Velocity Fluctuations . A gust 
has been shown to consist of a wide range of sizes of bursts or eddies 
of strong, high frequency fluctuations, separated by relatively quiet in¬ 
tervals. The cat’s paw itself is merely a burst on a lar^ time scale. 

The absolute magnitude of a single fluctuation is thus dependent not only 
on the instantaneous departure from adjacent values, but also on the de¬ 
parture of the fluctuation and its neighbors from similarly sized parcels, 
and on the departure of the mean magnitude of those patches from segments 
on yet larger scales; and so on up the range of scales to an averaging 
period covering the largest perturbations involved. Thus each scale of 
motion involves fluctuations whose magnitude depends both upon the 
strength of intel*nal instabilities, and upon the forcing of and interac¬ 
tions between the scale itself and larger scales or larger anomalies. 
Following Cramer (1955), let the instantanious magnitude of a random 
variate be ; and let the change of this variate caused by a for¬ 
cing be proportional both to the forcing and to the magnitude 

which had in turn been produced by the previous impulses. 

Then Si , c = 


and 




which in turn, in the limit of small change. 



For large N, the central limit theorem states that 


is normally 






61 . 


distributed. Hence has the lognormal distribution, 



- (JU(t-A') - 


') 


i X > X srd ^ 


where P ^ (x) is the probability '§ = x, CS"^ is the standard deviation 


of 3 , E(x) = ^ ^ ^ is appropriate normalizing or cut-off 

value of ^ 

The lognormal distribution can thus be expected to occur when a 
large number of impulses, or a large range of scales, interact to con¬ 
tribute to the increase of a variable at a rate proportional to its in- 
stantaieous value. For example, Rao and Narayanan (1971) have shown 
that the intervals between bursts in a turbulent boundary layer were log- 
normally distributed. This indicates that the scales of the bursts and 
intervals are determined by coupling with the scales of the boundary and 


of the outer or main flow as suggested by Mollo-Christensen (1970). 

Within a cat’s paw or gust, one variable which should act in a 
lognormal fashion is the turbulent energy dissipation, £ . Dissipation 

is primarily dependent upon the magnitude of the‘velocity fluctuation 


rate, which in turn depends strongly on the amount of ’rubbing’ or inter¬ 


action between different scale eddies. The dissipation is due to the 
working of molecular viscous stress against this coupled rate of strain 
or deformation of the fluid particles. Provided that the range of scales 
involved in the deformations is sufficiently large, say including inter¬ 
action between the localized high frequency instabilities, fluctuations 
on the scale of bursts and intervals, and similarly between groups of 
sub-bursts and bursts up to the scale of the full cat’s paw, dissipation 
should be approximately lognormally distributed. 

Another variable which should be lognormally distributed, and for 
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much the same reasons, is the magnitude of the squared time derivative 
of horizontal velocity. Time derivatives, like space derivatives, are en¬ 
hanced by the interactions between a variety of scales of eddies or in¬ 
stabilities. The interaction between high frequency velocity fluctuations 
and bursts has already been demonstrated in the preceeding section. Define 

the normalized squared magnitude of the velocity derivative as 

( 

> 

Chen (1971) has analyzed the distribution of with data taken 

from FLIP during BOMEX. He pointed out that squaring the derivative 
emphasizes the large fluctuations and minimizes the magnitude of small 
ones. (Gurvich and Yaglom [1967] showed that a lognormal distribution 
with large variance corresponds to a high probablity of values close to 
0 or much larger than the rms, but small probablities of intermediate 
values.) The natural logarithm of such a variable will have ^ 
as 'S O due to the ^ ’tail’ imposed by the rapid 

axis crossing of . Taking the average of a small number of 

y M/fc \ 

values (i.e., " X ^ ) wipes out this tail. In physical 

terms this corresponds to averaging the fluctuations over a very small 
but finite volume, or numerically to limiting the analysis to scales 
slightly larger than the cut-off imposed by the sampling rate. This aver¬ 
aging will increase the tendency toward lognormality, provided that 
the number of scales interacting is sufficiently large and that the averag¬ 
ing time is short enough to include the important structure. Chen’s 
data were from a roughly 4-minute period at mean heights of 3 and 31 meters 
above the sea surface. The resulting distribution was almost identically 
lognormal for averages over 8 data points. 
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Similar results from runs 2 and 6B of this study are shown in Figs. 
25 and 26. The burst/no burst/both character of the T.P. is noted in 
the figures. Fig. 25, with the data in histogram form (run 2, 2-second 
T.P.), clearly shows how averaging (over 5 points) removes the ’tail’ of 



the distribution. In Fig. 26, the cumulative probabilities of 


from both runs, 2-and 1-second T.P., are plotted on a normal probability 
scale. Purely lognormal distributions with comparable means and vari¬ 
ances are indicated by straight lines. 

Since only a relatively small number of points was used, the central 
limit theorem can not be strictly applied and none of the plots are 
strictly lognormal, though most of them show relatively minor departures. 
As shown best by the 1-second periods from the low stability run 2, per¬ 
iods of intense activity show a deficit of intermediate values, an excess 


of large and small ones. The quiet periods show the opposite effect, 
particularly in T.P. 6, taken after the passage of the gust (past the 
end of Fig. 2). In this very quiet T.P. the fluctuations are more nearly 
normally than lognormally distributed. The 2-second T.P. show these 
effects less clearly due to averaging over both bursts and intervals, 
but the tendency toward normality is clear in the quiescent T.P. 5. The 
run 6B data (high bulk stability) shows similar results, with intermediate 
deficits during bursts, excesses during intervals. One major exception 
is T.P. 2, from the leading edge of the cat’s-paw. It was shown in the 
preceeding section that in this time period the major effect of the 
stress fluctuations was the local tumbling and destabilization of the 
boundary layer. The excess of intermediate departures in this distribu¬ 
tion is further evidence of the initial effectiveness of the buoyancy 
forces in retarding momentum transfer and limiting fluctuation strength. 
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The distribution of the squared velocity derivative is thus an ef¬ 
fective indicator of the range of scales involved in boundary layer tur¬ 
bulence. Interactions are strongest and involve the largest number of 
scales during the bursts. During quiescent intervals and particularly 
between gusts, fewer scales interact; coupling is weak or non-existant, 
and tends to have a normal distribution, indicative of a lack of 

interaction. 

It should be noted that for isotropic, homogeneous turbulence,i 
is related directly to dissipation by use of Taylor's frozen turbulence 
hypothesis, ^ r - “U 3^ . Chen in fact regarded his data in 

these terms. Ruggles (1969) has demonstrated that the frozen turbulence 
hypothesis is quantitatively inapplicable in a case even less restrictive 
than this study; and indeed the above results indicate that it is not 

I 

necessary to so restrict the assumptions on the nature of the flow in or¬ 
der to account for the lognormal behavior of fluctuation magnitude. The 
well-known intermittency of the small scale structure of turbulence, the 
consequent large dispersion of £ , and the dependence of the probability 

structure of 6 upon larger scales (Obukhov [1962], Landau [1959]) 
have been theoretically applied to the statistics of turbulence only for 
the locally isotropic, homogeneous case (Kolmogorov [1962]). The present 
development has avoided these restrictions and emphasized strong non¬ 
linear coupling between scales as the causative phenomenon. Further, 
it has been demonstrated that local departures from lognormality are 
also indicative of the relative strength of coupling of the turbulent 
components, and that many features of the boundary layer and transfer 
processes, such as the spacing of bursts and intervals and the frequency 
and strength of local velocity fluctuations, are controlled by this 
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coupling. 


B.2. An Example of Weak Resonant Wave Interactions . 


Nonlinear wave-wave interactions have been theoretically studied in con¬ 
siderable detail during the last 10 years (Kinsman [1965] provides a 
good summary of the major concepts). Energy exchange among gravity 
waves of different wave numbers is only seen when the perturbation ana¬ 
lysis extends to third order; it involves triads of primary waves inter¬ 
acting to feed energy into a fouth component. For waves in the cap¬ 
illary-gravity range, the analysis shows that interaction occurs at 
second order (as well as at third order); two waves interact to feed 
energy into a third (McGoldrick [1965], Simmons [1969]). This weak 
resonant interaction has been demonstrated in the laboratory (McGoldrick 
[1970]). 

McGoldrick (1965) has discussed the implications of this energy 
transfer with respect to wind-wave generation. Of particular importance 
is the fact that the interaction times involved in the growth of the re¬ 
sonant wave are only 2-3 periods of the initial waves creating it. Be¬ 
cause of this short interaction tim^, the amount of energy that the wave 
field can receive from the wind at a given frequency is determined not 
just by breaking but also by nonlinear transfer to other components. 

Of particular interest is the * self-resonant’ wave, of wave number 



( T = surface tension/density). This is commonly 


known as Wilton's wave; he noted the presence of a singularity at this 
wave number in his (1915) perturbation analysis. The general require¬ 
ments for resonance at 2nd order are k 2 + k^ = 0, ± W 2 + = 0, 


and ^ ^ i 


where w 
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fulfills this relationship for = -1/2 k 2 , or ' 




Simmons noted that waves shorter than k 2 are hydrodynamically unstable 
at 2nd order; further, k 2 is the longest wave which can interact with 
two longer waves. Waves at k^ are extremely unstable and can occur only 
as a transient state; they can not continuously exist as free surface 
waves and once generated will rapidly transfer their energy to the k 2 

^2 


I—I . 3 2 

waves. For clean water, i. = 75 cm /sec and k^ and k^ roughly 


correspond to frequencies of 10 and 20 hz, respectively. 

Fig. 27 shows the wave autospectra from the 2-second segment ana¬ 
lysis of run 6B.(80% confidence limits are roughly 1 decade wide). All 
5 T.P. display the peak of the major wave field at around 1 hz. T.P. 1, 
from before the gust, is rather smooth in the band above 3 hz. T.P. 2, 
near the leading^edge, shows a slight deficit around 9.5 hz, followed by 
a levelling off to 14 hz, a rapid drop to 19 hz, and a slight increase 
around 20 hz. During the ’active' T.P. 3 and the major burst in T.P. 4, 
spectral slope is less steep and the spectrum shows a slight bulge near 
10 hz (i.e., near the frequency where leveled off and reversed 

slope); in T.P. 4, there is also a Local maximum just above 20 hz. In 
T.P. 5, immediately after the burst, there is a significant deficit 
near 10 hz, and again a local minimum near 19 and. maximum just above 
20 hz. This sequence of events, particularly the change between T.P. 4 
and 5, is strongly suggestive of the picture given by McGoldrick. Dur¬ 
ing generation the k^ waves are not only present but near saturation 
because of the energy input from the wind; as soon as the forcing dis¬ 
appears, so do the extremely unstable k^ waves. Their energy is 'drained' 
to higher frequencies. Fig. 28 shows a similar example from run 2, 
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where the mean wave field was significantly larger and the results con¬ 
sequently less clear. Here the major burst was in T.P. 3, with lesser 
activity in T.P. 4 and quiescent in T.P. 5.’ In T.P. 5 the spectral 
deficits occur at 10 and 18 hz, with the local maximum around 19.5 hz. 

The slight frequency shifts of the postulated and waves are 
readily attributable to contamination of the Charles River. The very 
fact that the predicted effects of weak resonant interactions stand out 
at all in other than a ’carefully tuned' laboratory experiment, is 
strong evidence of their importance in wave spectral development. Ev¬ 
idence that perturbation analysis is suitable for predicting continuum 
nonlinear effects in major wave fields was recently provided by the 
bispectral analyses of Garrett (1970). Nonlinear transfer of energy 

is thus apparently of importance at all stages of spectral development. 

c*. 

C. A Further Note on Time Averaging . The results of the preceed- 
ing sections have emphasized the importance of high frequency fluctuations, 
the intermittent nature of generation, and the coupling between scales. 
Averages over time scales comparable to the entire gust or to a combina¬ 
tion of burst and interval have been shown to obscure the selective 
high growth rates associated with the local instabilities, the effec¬ 
tiveness of wavelets in coupling the wind and the waves, and the nature 
of interactions between scales. One further example will serve to il¬ 
lustrate the critical importance of proper averaging in* the analysis of 
wind-wave generation data. 

Consider a typical cat's paw from run 2. Typical mean velocity 
values for different scales can be very roughly taken as 








72. 


u 

Z noors 



T'gus't ■ *“5 4 

boryt 

2jO 

> 



Similarly, temperature means may be estimated, with 
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Next, assuming rms turbulence levels of 10% (a very low estimate for 
a gust or burst), 3.75° C 4.25°, 18 m/sec ^*burst^ m/sec. 

As demonstrated in Section A, the important fluctuations in heat 
and momentum transfer and wave and turbulence generation are = w*T* 

and T - -y:? u*w’ (assuming symmetry in v’ and V = 0) . Then 

U=<U>+u’, W=<^W>+w’, T=^T>+T* 

and since 

W > and^wSare taken = 0 over whatever averaging period is chosen, 

u’w* =w*(U-^U>) ‘ and w’T’ =w'(T-CT>). 

Because of the impressed zero mean of W, the long term mean fluxes 
will be independent of the choice of U > and <^T> ; however, both the 
magnitude and the sense (sign) of the instantaneous fluxes (or of very 
short averages) are strongly affected by this choice. For the example 
above, choosing U = <U> ^ ” ^^^gust’ plotting a joint histo¬ 
gram of the resultant values of ZT and for the burst, the result 

will appear as in Fig. 29B (Fig. 29 is from run 2, 2-second T.P. 2); while 
if the means are chosen as the burst means, results are as shown in Fig. 29A. 
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When viewed against the general gust background it appears as if the ^ 
and fluctuations are negatively correlated (of opposing effect in 

turbulence production) and very strong; when viewed against the 2-second 
background, mostly burst, the fluctutations are relatively much smaller 
and show a weak positive correlation. Similarly, when the results are 
applied to the turbulent kinetic energy balance, the effects on production, 
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and hence on the stability structure of the boundary layer, are opposite. 


This is an extrememly strai^tforward but seldom applied result. The 
proper choice of averaging period must be made on the basis of the im¬ 
portant scales in the process under consideration. For wave generation 


and air-sea momentum transfer, these have been shown to be the scales of 
and within the bursts and the main wave field. Consideration of the wave 
generation process in terms of long term mean effects can totally obscure 
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the physical realities. 


Summary of Results and Conclusions 

The primary facts to emerge from this study may be summarized as fol¬ 
lows , 

1. Cat's paws are highly intermittent phenomena that them¬ 
selves are composed of localized bursts separated by quiescent intervals. 

The bursts contain high frequency fluctuations of stress and buoyancy 
3-4 orders of magnitude larger than that observed in the mean. The var¬ 
iability in time was shown in Fig. 2 - Fig. 7; representative magnitudes 
are shown there and in Fig. 8. 

2. The bursts are associated with the generation of high 
frequency waves and turbulence. This is best demonstrated by the se¬ 
quence of events shown in Fig. 23. Participation of the fluctuating 
boundary enhances the magnitude of the fluctuations, as shown by a com¬ 
parison of the results of this study with measurements over land. 

3. The high frequency instabilities not only have large growth 
rates on their internal scales, but also are associated with the growth of 
lower frequncy components. Both high and low frequency fluctuations are 
strongly enhanced in the bursts, as shown by the changes in magnitude and 
shape of autospectra (Fig. 9) and on the detailed SC-4020 plots. This is 
not only true for couplings within but also for couplings between media as 
demonstrated by the diffeienoes in coherence between stress and surface 
waves in the 2- and 1-second time periods (especially T.P. 4, Figs. 21 and 
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24). The bursts are also associated with the development of’ shear layers 
in the wind profile and with large values of turbulence production (see 
Fig. 10), a further indication of the importance of localized instabilities. 

4. The boundary layer is mixed and tunbled by the instabil¬ 
ity fluctuations in the bursts. The effect of strong temperature strati¬ 
fication in the boundary layer is to delay the onset of wave generation by 
requiring intial utilization of available turbulent kinetic energy for 
destabilization. Buoyancy forces are very effective in degrading the 
magnitude and effectiveness of momentum flux, and quickly cause restabil¬ 
ization after passage of the gust. 

The general mixing and boundary layer inversion patterns are shown 
in the joint histograms of U - T’ (Fig. 11) and (Figs. 13, 15) and 

in the phase reversals in the coherence of T’ - ^ and “c -jS (Figs. 12, 14). 
The mixing associated with the bursts, and the partial restabilization 
during intervals, is shown by the changes in correlation between ~o and 
/3 in Fig. 20. The delay of generation attributable to the utilization 
of kinetic energy for mixing is demonstrated in Figs. 22-24, where the 
burst at the leading edge (T.P. 2), uncoupled to the waves, has little 
effect upon the water surface but effectively inverts the boundary layer. 

The effects of stability have also been demonstrated by the ex¬ 
cess of intermediate sized velocity fluctuations in the leading edge 
(Fig. 26, runs T.P. 2), by the rapid rise of T after gust passage (Figs. 

11, 15), and by the temperature increases and decreases associated with 
T and y6 fluctuations of different relative strength (e.g.. Fig. 6). 

5. The coupling between scales is strongly nonlinear, as 
illustrated by the lognormality of the distribution of squared velocity 
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fluctuations during bursts (Figs. 25, 26), the strong association of high 
frequency instabilities with large magnitude velocity, stress, and buoy¬ 
ancy perturbations (Figs. 2-7 and 9), and the association of burst length 
with the period of the major underlying wave components (Figs. 2, 4, 6). 

6. Weakly resonant nonlinear interactions are effective in 
transferring energy among capillary-gravity waves. This was demonstrated 
by the change in the and k^ components (Figs. 27, 28) between*’genera¬ 
tion’ (burst) and ’detritus’ conditions. 

7. The important momentum transfer processes are of short 

duration, highly nonlinear, frequency selective, and strongly associated 
with intermittent high frequency instabilities and coupling between scales. 
This was demonstrated by the variations in coherence patterns between 

different averaging times and burst/no burst conditions (Figs. 21, 24), by 
the simultaneous occurrence of strong Reynolds stress and capillary wave 
activity,(Fig. 23), and by the selectivity of coherence and spectral change 
(Figs. 24, 27, 28). 

8. The differences between ’generation’ and ’detritus’ can be 

described by the changes in coherence between variables and in the proba¬ 
bility distribution of velocity fluctuations. Detritus coherence is char¬ 
acterized by low magnitude (^ ) and/or by continuous phase ( © ) slippage 
in any coherent bands; during strong coupling and generation, T is high 
and & is nearly constant in the frequency bands that participate most 
strongly in the transfers (Figs. 19, 21, 24). Velocity fluctuations 
(measured as ^ ) in detritus have excesses of intermediate values and 

approach normal distributions; those from generation periods are more 
nearly lognormally distributed with excesses of very large values and 
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values near zero. These differences are due to the enhanced' range of scales 
and strength of the couplings involved in wave generation and turbulent 
interactions. 

2 

The intermittency switch function (X = 1, ^ 100 dynes/cm , I = 0, 

2 

f < 100 dynes/cm ), shown only in Fig. 5 but computed in all runs and 
used for selection of the 2 and 1 second time periods, appears to be a 
suitable first approximation for a conditional sampling switch. When 
smoothed over a time period ^0.5 seconds it provided a reliable indicator 
of the presence of a burst and instability-associated generation. 

9. There is a rough, qualitative similarity between scales 
of motion (e.g., gust and burst), events of different magnitude (Fig. 7), 
and between turbulence patterns over land and water surfaces. The similar¬ 
ity lies in the presence of a large range of scales. The significant dif¬ 
ferences are caused by the strength of coupling between these scales; the 
larger the number of participating scales and media, the stronger the coupl¬ 
ing and the more efficient the associated heat and momentum fluxes. These 
features were demonstrated by a comparison of flux measurements over land 
and water, and illustrated in the time plots of the variables. The simul¬ 
taneous enhancement of frequency range and fluctuation magnitude in 
bursts is perhaps best displayed in Fig. 7, where the ’instantaneous’ 
magnitudes are based on the 10 second means. 

It was also demonstrated (Section C, Fig. 29) that the effect of t* 
and ^ fluctuations can be different when their magnitude and sense are 
evaluated over different scales; e.g., a burst may contribute to boundary 
layer instability as seen from the long term mean, while the local and 
immediate effect of the fluctuations is to temporarily reestablish a 
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portion of the initial temperature stratification. Averaging’ over time 
scales other than those of importance in the momentum transfer process 
can lead to misinterpretation. 

10. Evaluation of generation in terms of the ^mean’ can be 
misleading and obscures the physics of the transfers. Since much wave and 
turbulence generation does not occur on the mean but in .localized, fre¬ 
quency selective, high growth rate instabilities (as depicted in Figs. 23, 
24, and 27), and since the mean itself is a seldom if ever assumed state 
(as shown in the time series diagrams [Figs. 2-7] and the histograms 
[Fig. 8]), representation of the generation mechanisms only in these terms 
is unsatisfactory. 

On the basis of these results it is concluded that microscales play 
a predominant role in wave generation. High frequency waves and turbulence 
are intermittently generated during high growth rate, local instabilities, 
and the wavelets subsequently act as drag elements to couple the longer waves 
to the wind. Further, the high frequency instabilities aid the growth of 
similar, lower frequency fluctuations. This combination of events leads 
to spectral growth rates much larger than those predicted on the basis of 
linear and weakly tionlinear theories. 
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III. THE ROLE OF MICROSTRUCTURES IN WIND^WAVE SPECTRAL DEVELOPMENT 


Abstract 


Measurements of the wind, wave and surface current fields were made 
from a 38.5 meter spar buoy at the mouth of a tidal channel. Statistical 
analysis of velocity profiles showed that both air and water boundary lay¬ 
ers were predominantly logarithmic (neutral stability),.and confirmed 
Ruggles^ (1970) finding of anomalous behavior of atmospheric profile para¬ 
meters at critical wind speeds. A lack of simple relationships among pro¬ 
file measures indicates that wind-wave coupling can be only very poorly 
characterized on the basis of means, that the time history of spectral 
evolution is crucial to an understanding of the physics governing the mo¬ 
mentum transfer, and that more than one wave generation regime exists. 

A two-day history of spectral development, which did not include the 
major profile anomalies, was analyzed to determine the nature and scales 
of momentum transfer processes. Wave generation was found to occur pre¬ 
dominantly in short, intermittent, burst-type events, and to involve 
strpng nonlinear coupling of wind and wave fields with associated rapid 
growth of the high frequency components. These important events involve 
coupling between wave groups in wind and surface wave fields and strong 
interactions between the media and between frequency bands in both fields. 
Tidal currents modulate the wave field and cause imbalance of wind-wave 
equilibrium conditions. 


Introduction 

Microscale studies in cat’s paws have demonstrated the importance of 
intermittent, strong interactions between scales of motion and between at¬ 
mospheric turbulence and the sea surface in the generation of high frequen¬ 
cy waves. Similar types of interactions are also significant in the de¬ 
velopment of the major portion of the wind-wave spectrum. 

g 

In the open ocean, typical Reynolds numbers are on the order of 10 ; 
local wave generation occurs under the influence of surface drift, tidal 
currents, propagating swell, and stability, precipitation, wind speed 
and direction shifts associated with large weather systems. The large 
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number of participating scales and phenomena can all interact, with both 
strong and weak coupling, to have a significant effect on spectral de¬ 
velopment. It will be shown that there are separate regimes of wave gen¬ 
eration, and that the time history of events is of crucial importance in 
development of the wave field. 

Most previous open ocean studies of spectral development have dealt 
with mean or long time average wave energy growth and have examined the 
influence of only atmospheric conditions upon the wave field. Airborne 
measurements (e.g., Barnett and Wilkerson, 1967; Schule et al, 1971) and 
other laboratory and oceanic studies which extended over a range of fetch 
(e.g., Snyder and Cox, 1966; Sutherland, 1967; Mitsuyasu, 1969) have 
dealt extensively with the major features of spectral development with 
little or no experimental reference to the causative momentum transfer 
mechanisms. From these types of studies spectral growth rates, the fre¬ 
quency shift of the spectral peak, the existence of an equilibrium range 
for high frequencies, and the overshoot/undershoot phenomena have become 
well known features. 

Studies from stationary buoys or platforms (e.g., Seesholtz, 1968; 
Ruggles, 1969) have evaluated growth rates in terms of wind velocity pro¬ 
file characteristics with minimal consideration of the effects of the 
underlying current and subsurface boundary layer conditions. Such studies 
have shown that the atmospheric boundary layer is essentially logarithmic 
(at least for near-neutral stability) and that spectral growth rates are 
substantially larger than those predicted by the inviscid laminar Miles 
theory (Miles, 1957). 

Investigations of the direct interplay between wind and waves have 
concentrated on the transfer mechanisms suggested by this theory. One of 
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its critical elements is that turbulent and wave induced Reynolds stresses 
do not interact, and that momentum transfer occurs by means of an insta¬ 
bility in the mean shear flow. The momentum is extracted from the mean 
flow in the critical layer and transmitted downwards by wave induced Rey¬ 
nolds stresses. 

The recent work of Kendall (1970) demonstrated that a fixed wave 
pattern caused significant modulation of the turbulent structure^ but 
that the induced fluctuations did not have a constant phase relation to 
the waves. Davis (1970) and Stewart (1970) theoretically and experimen¬ 
tally investigated the vertical structure of Reynolds stress.and found 
poor quantitative agreement with the wave-turbulence Reynolds stress pro¬ 
file predicted by instability theories. Wave induced Reynolds stresses 
accounted for only about half of the momentum carried to the wave surface 
by the pressure component in phase v/ith wave slope. 

Ruggles (1969), Davidson (1970), and Shemdin and Lai (1970) have all 
demonstrated the existence of wave induced perturbations in the wind 
field and have shown that these wave induced Reynolds stresses do interact 
with turbulent stresses. There is thus substantial evidence that inter¬ 
action between atmospheric turbulence and surface waves is one of the cru¬ 
cial mechanisms involved in the transfer of momentum to the waves. This 
problem is being studied by Saeger and Reynolds (1971); their initial re¬ 
sults indicate that the inviscid model, by neglecting the important struc¬ 
ture of the critical layer, gives very poor predictions of the Reynolds 
stress and hence of momentum flux to the waves.. 

The cat’s paw studies have demonstrated the importance of detailed, 
near-surface measurements in determining the interactions among forces. 
Dobson (1969,1971) recently made the first successful detailed field 
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measurements of atmospheric pressure on wind generated waves. His im¬ 
portant results showed that the observed phase shifts were an order of 
magnitude larger than predicted by Miles' theory, and that energy and 
momentum flux spectra were strongly peaked at- or above the peak of the 
locally generated waves. Further, the inviscid laminar Miles model 
underpredicted increases in wave energy by factors of 5-8. These facts 
suggested a strong interaction between wave induced and turbulent Rey¬ 
nolds stresses. Dobson also noted that energy flux and wave supported 
momentum flux were negligible above relatively low frequencies, contrary 
to an assumption that much of the momentum transferred goes into short 
waves (Wu, 1969; and the cat's paw and present studies) but not contra¬ 
dicting the suggestion that short waves aid the growth of lower frequency 
instabilities by nonlinear coupling between scales of motion. Near sur¬ 
face pressure measurements over waves made by Elliott (1970) confirmed 
Dobson’s pressure shift findings. Elliott noted that pressure lagged 
the wave by about 135° during generating> and by 180° during non-genera¬ 
ting conditions. There was no vertical phase shift at heights above the 
crests. Elliott's work also suggested strong interactions between wave 
induced components and normal boundary layer turbulence. 

The basic conclusion from detailed investigations of near-naviface 
structure and air-sea coupling is that turbulence and waves interact 
strongly and intermittently. Only a relatively small portion of the 
wave growth is due to processes which act on the 'mean', and these pro¬ 
cesses are themselves strongly dependent upon the presence and location 
of high frequency components and associated strong nonlinear coupling. 

The structure of the underwater boundary layer is also of major 
importance in transfer processes. It has been repeatedly demonstrated 
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that the magnitudes as well as the mechanisms involved in transfer of mo-- 
mentum, heat, and moisture across an interface are radically different for 
flat plates or solid surfaces and fluctuating or moving surfaces (e.g., 
Chambers et al, 1970; Mangaralla et al, 1970). Over water, the processes 
are different for the cases of wind waves alone and wind waves with the 
added presence of swell and surface drift. 

Few experimental studies have considered the effect of the water 
near-surface layer motions upon wind-wave coupling, although the weak in¬ 
teractions between waves and currents or very long swell are well known 
theoretically (e.g., Longuet-Higgins and Stewart, 1960, 1961). The 
available experimental evidence on the structure of and interactions in 
the water boundary layer is somewhat contradictory. Shonting (1967, 

1968) noted that u-w phase relations were close to those predicted by 
classical theory, with no phase shift with depth and only minor interac¬ 
tion between wind waves and swell. Shonting, Stewart and Grant (1962) 
and Yefimov and Khristoforov (1969) have all found that wave stress is 
strongly concentrated very near the surface, and that essentially all 
dissipation occurs above the trough. The transfer between high and 
low frequencies caused by dissipation (breaking) and viscous decay of 
short waves (Longuet-Higgins, 1969b; Hasselman, 1971) is also a surface 
phenomenon, and involves non-linear coupling between scales. Shemdin 
and Lai (1970) found that the wind-induced drift in a wind/water tunnel 
had a logarithmic profile, and that the surface stresses computed from 
the air and water profile parameters were roughly equivalent. Intro¬ 
duction of mechanical waves disrupted the wind driven current. The 
couplings between the adjacent boundary layers and the interface are thus 
as yet imperfectly described. 
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This experiment was designed to investigate the nature and time 
scales of air-sea coupling during wave generation. It was conducted 
during a variety of wind, weather and current conditions. The analysis 
was divided into three segments. First, a study was made of the air and 
water boundary layers in terms of a logarithmic profile model. This 
section extends and complements the investigations of Ruggles (1970) 
and Groscup (1971), and describes the general environmental background 
of the more detailed analyses. Second, the time history of spectral 
growth during a two day period of increasing winds was investigated. In 
this section the general nature of coupling between wind, wave and current 
fields is described. Finally, segments of the two day record were studied 
in detail to emphasize the importance of strong, intermittent generation 
of high frequency waves and resultant wave-turbulence couplings that were 
demonstrated in the cat^s paw study. Generation ^events’ are described 
and their influence on the main wave field is displayed. 

The major objectives of the study were to demonstrate the wide 
range of interactions involved in wave generation and spectral develop¬ 
ment, and to emphasize the importance of high frequency ^microscales^ 
and intermittency in the momentum transfer process. 


The Experiment 

Measurements of wind and water speed profiles and wave height were 
taken from a stable, 38.5m spar buoy (Mollo-Christensen and Dorman, 1971) 
(Fig. 1). The buoy was slack moored in 31.1m of water (MLLW) and faced 
into the tidal current (NE flood, SW ebb, 1.5m/sec maximum current). 

The experiment ran from July through October, 1970. The site was at 
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igure 1. M.I.T. Air-Sea Interaction Buoy 
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the mouth of a tidal channel (Vineyard Sound, Massachusetts); minimum 
fetch to the S and SW (direction of prevailing winds and sea breeze) was 
50km. 

Instruments included 5 cup anemometers at heights of 2.4m to 11.4m, 

3 ducted-propeller (5cm diameter) current meters at depths of -1.5 to 
-4.7m, and two capacitance type wave gauges with 4.6m sensor length. The 
primary wave gauge was mounted directly over the current meters, 2.5m 
cross-stream from the anemometers. Thermistors were attached to arms of 
the velocity transducers at +5 and -2.5m. Buoy tilt was sensed with a 
pendulum pot. 

Signals from the transducers were mixed, telemetered to a shore sup¬ 
port facility on Cuttyhunk Island, Massachusetts and recorded in analog 
mode. Time, voice log, and 25kc clock signals were transmitted and re¬ 
corded for use in synchronization. Velocity information was in pulse 
format, wave gauge signals a variable frequency, and thermistor data a 
variable D.C. voltage. Data processing was done at M.I.T. and used a 
combination of analog and digital methods. Digitization rate for the 
time history analysis was lOhz; velocities were pre-filtered at 2.5hz, 
wave gauge signals at 5hz, to avoid aliasing. 

A total of 146, 10 minute ’mean condition’ runs were made. There 
were 5 series of runs, with data taken once each hour for periods ranging 
from 6 to 53 hours. 109 of the runs contained data suitable for analysis. 
Winds and waves were from the S-SW during almost all runs, with 10 min¬ 
ute mean wind speeds (at 10m) of 1.5 to 12.5m/sec. Because of the large 
fetch the atmospheric boundary layer was nearly neutrally stable during 
all runs (AT air-sea 3®C and -.02 < Ri(bulk^ ^ .05). Similarly, the 
water surface layers were generally well mixed down to a depth of at 
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least 8-lOin. 

The analysis was conducted in three parts. First, a distribution 
analysis was made of the parameters obtained from a comparison of velocity 
data with a neutral stability logarithmic profile model. Of the 109, 

10 minute anemometer runs, 96 (88%) had correlations 0.94 with the 
model. Of the 51 minimeter profiles used (from the longest nearly con¬ 
tinuous segment, runs 109-159, also used in the time history analysis), 

48 or 94% were logarithmic to within the required contraint; 49 of the 
anemometer profiles from the same period were logarithmic. Much of the 
analysis was based on 30 second averaging times. For the anemometers 
there were 2467 samples (81% logarithmic; 1057 samples with 91% loga¬ 
rithmic for runs 109-159), for the minimeters 1060 samples (85% loga¬ 
rithmic). Ruggles ( 1969, 1970) has discussed the significance of the 
correlation coefficient; Groscup (1971) has described the increasing 
tendency to a logarithmic profile with increasing averaging time for the 
runs discussed herein. Only the runs satisfying the ^ 0.94 criterion 
were used in the histogram analysis, since the parameters are computed 
from the assumed model. This procedure is really only satisfactory for 
parameterization and description of the general background conditions 
and generation regimes, since as shown in the cat’s paw study strong 
intermittent events are associated with large profile deflections and the 
development of unstable shear layers. The approach to a logarithmic pro¬ 
file with increasing averaging time is merely another indication that long 
term averages can obscure the physics of the transfer processes. 

The longest set of runs was then analyzed using standard time series 
techniques. The data is presented in terms of spectral growth and the 
related changes with time of atmospheric and water variables, and the 
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coherence-phase relationship between horizontal wind speed and waves. 
Spectral estimates have 6 to 20 degrees of freedom Cdof); the 80% and 
90% confidence limits of the distribution (at $ - 10^, n integer) 

are shown on the spectral growth and spectral' density estimate diagrams. 
Garrett (1970) has recently demonstrated the independence and normal dis¬ 
tribution of the real and imaginary parts of the Fourier coefficients im¬ 
plicit in the use of this measure for mean wave fields. The confidence 
limits are applicable only to each individual realization of the spectrum, 
and do not provide estimates of variance of the wavefield .for times longer 
than the overall averaging period; e.g., the 30 second spectra are strong¬ 
ly affected by the presence or absence of swell and major wave groups. 
Coherence magnitude ( ) was tested on the assumption that the variates 

were random with a bivariate gaussian distribution (Amos and Koopmans, 
1963). (As demonstrated in the cat’s paw study and further displayed 
herein, this poor assumption for strongly coupled variables; un¬ 

fortunately it is the only readily available method). The results may be 
generally stated (for 5 minute time segments) as, at the .05 level of 
significance, the probability that the hypothesis = 0 is incorrect 
for a measured ^ ^ 0.5 is 97%. (For the 30 second segments, this 

probability has dropped to 70%). 

Finally, several 10 minute segments from the start, middle, and end 
of the set of 51 runs were analyzed using 30 second averaging periods. 

The data is presented in a similar fashion, with the added utilization 
of analog strip chart displays of ’events’. 
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Results and Discussion 


A. Boundary Layer Profiles . With the assumption that the averaged 
stress ( T ) in the near surface region is constant with height, a charac¬ 
teristic velocity of the turbulent boundary layer may be defined, 

r J-U V’ . 


This is termed the friction velocity. Using similarity arguments, for a 
neutrally stable atmosphere the only characteristic length is height. 

The resulting profile equation is 


u(0 = 


where ^ 


is Von Karman’s constant - 0.40 and the roughness length, is general¬ 
ly not added to the height in the numerator (see, e.g., Lumley and Panofsky, 


1964). 


There is ^neral agreement that this simple model gives a good 

description of the mean velocity profile in conditions of neutral stability 

over both land and water, but there is still considerable debate on the 

meaning of and relationships between the profile parameters (U. and z ) 

« o 

and other variables computed from the model. From the development of the 
model, and z^ are independent similarity parameters characterizing 
the boundary layer; however, Sadeh, et al, (1971) have demonstrated that 
once the roughness elements are sufficiently high compared to boundary 
layer thickness, U. and z become local variables. 

Measurements of the profile parameters over waves have been re¬ 
viewed and discussed in detail by Wu (1969a, 1969b, 1970) and Ruggles 
(1969, 1970). Wu states that the wind stress (7* ) is supported by 

form drag provided by wavelets superimposed on the larger waves; air flow 
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separates from the wavelets whose phase speed is less than U^. becomes 

a measure of the height of these roughness elements (roughness height 
= 30 z^) and is related to friction velocity by Charnock’s relationship, 
zj ( ) = constant. This is manifested by an increase in the 

drag coefficient (Cd^ = = ( ^*'/ ) ) with U up to a ’sat¬ 

uration* value, U = 15 m/sec. Above 15 m/sec the wind speed is larger 
than the average wave phase speed and Cd and z^ become nearly constant. 
Effectively, the Charnock relation implies that a single parameter can 
be used to describe the flow and thus the momentum balance of waves and 
atmosphere. The development of the profile model on the other hand as¬ 
sumed that and z^ are statistically independent. 

Ruggles showed that over waves, the wind speed is more appropriate¬ 
ly represented as the sum of a mean profile component (supporting atmos¬ 
pheric dynamics), and a term representing the correlation between and 
z^ (from wind-wave interaction). The ’correlation’ wind component repre¬ 
sents that port(i^n of atmospheric momentum required to support air-sea 
interaction dynamics. It is roughly 10% to 20% of above 4 m/sec, 

and negative below 2 m/sec (indicating a minimum wind speed for generation). 
Ruggles also found marked changes in z^ and Cd, with excessively high peaks 
and excessively large variances, at around = 2, 4, and 8 m/sec. These 

anomalies correspond to the division points of the Beaufort scale, i.e., 
to marked changes in the physical appearance of the sea surface. Similar 
anomalies are shown by Wu, who attributes the changes to the differences 
of flow in aerodynamically smooth, transition, and rough regimes. Hicks 
and Dyer (1970) found a simpler relationship, with Cd = .0011 and aero¬ 
dynamically smooth flow below 6 m/sec, and Gd = .0015 with rough flow for 
Ufo = 10 m/sec. They did not report any anomalies. 
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Additional evidence of anomalous wind speeds was given in the work 
of Groscup (1971) and Mollo--Christensen and Groscup (1971) (based on the 
data of this study). They further demonstrated a bi- and tri-modal dis¬ 
tribution of in the vicinity of the = 8 m/sec discontinuity, and 

a bimodal distribution of U. at or near the z and Cd discontinuities at 

* o 

4 and 8 m/sec. The anomalies are indicative of enhanced air-sea inter¬ 
action at specific wind speeds and of changes in the nature of wind-wave 
coupling. The large variances and bimodal distributions are evidence of 
enhanced nonlinear coupling at critical wind speed/surface‘condition values. 

Fig. 2 shows the distribution (smoothed histograms) of 30 second 
average wind speed and of atmospheric logarithmic profile parameters from 
this study, both for all runs and for the subset 109-159 (those used in the 
time history analysis). Fig. 3 through Fig. 6 show selected joint histo¬ 
grams of the variables. 

^10 found from the profile model, using the computed values of 
and z . The drag coefficient, Cd, was determined for both the computed 

3 

10 m and measured 11.4 m speeds. is a measure of turbulent kinetic 

energy production, 

— -- - 

and Ou is the variance of velocity in the profile, 

The histograms show that a large percentage of the runs had very 
small roughness length ( (^o')^ - 7.^ or ^ y values of 

-3 

z^ <1 10 cm were lumped into the In [z^] = -10 category so that sta¬ 

tistics were not computed). Stress values are in rough agreement with 
the long term averages found by other investigators; the mean is largest 
for the continuous strong winds of runs 109 - 159, and decreases with 
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increasing averaging time. The values of Cd are comparable for the 
measured and computed heights; Ruggles (1969) computed Cd^ for 5 anemometer 
heights and 10 m, and found similar results> with as expected a slight 
decrease with increasing z. The Cd values are concentrated in the normal 
range for runs 109 - 159, but when all runs are considered together, the 
30-second and 10-minute means are excessively large and there are numerous 
exceptionally high drag values. These major, average changes in the drag 
effect of waves on the wind suggest the existence of more than one wave 
generation regime. , Cci,o :> and (and thus and r ) all have 

decreasing mean values with increasing averaging times, indicating that 
profile distortion and resulting turbulence and wave generation occur in¬ 
frequently; the short, strong, important anomalies are obscured by long 
term averaging. 

Fig. 3 shows JH (U^Q, U^), JH(U^q, In(z^) ), and JH(U^q, Cd^^) for the 

30-second average profiles. In this figure, anomalies of U. ln(z ), and 

Cd^Q all occur at roughly 2.5, 3.5, 6.4, 7.7, 10, 11, and 12.2 m/sec (they 

are particularly striking at 2.5, 3.5, and 7.7 m/sec, i.e., near those 

foxmd by Ruggles). The anomalies are in the form of exceptionally large 

variances and wide departures of the interval means from any postulated 

simple fit of the data. The values are in general quite scattered and 

2 

poorly represented by means. Also shown is JH(U^ /g, K) where K = 30 
is height of the roughness elements (wavelets).as postulated by Wu. It 
is obvious that a simple relation such as suggested by Charnok and found 
by Wu does not hold here. There are bimodal distributions and anomalies 
in both directions. A single parameter can not characterize the profile, 
because of the interplay between atmospheric dynamics and the wave field. 
None of the data can be fit with a linear relationship; in particular. 
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Figure 3. 
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Ruggles*s suggested = 0.04 U^q intersects the minima of the mean of 


10^ 

Fig. 4 shows ^^10^ computed using the formulae of Wu (1969b) 


for 'breeze* and 'light wind* regimes, 
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The anomalies at 2.5 - 3.5, 7.5, and 12 m/sec again appear, but much less 
strongly than in direct computation from the profile model; in addition, 
the mean value of Cd computed using this technique is excessively small. 
Large scatter was found in joint histograms of other relationships sugges¬ 
ted by Wu. The failures of these postulated simple models on the relative¬ 
ly short (30 sec) averaging times indicates that the dynamics of wind flow 
over open ocean wave fields at large Reynolds numbers can be only poorly 
described in terms of long time averages, even when the significance of 
wavelets is considered. Wu's contention that the small surface components 
support the majority of wind stress is not disproven by these data; the 
results merely indicate that the dynamics of the wind-wave coupling is 
considerably more complex than that caused by just form drag and air flow 
separation. 

Fig. 4 also displays 30 second JH( r , In(z^) ), JH(Cd^Q, In(z^) ) and 

JH( O"o /k? )• The T -ln(z^) relationship shows a distinct biraodality, 

with generally increasing stress with roughness length and stress maxima 

for ^ 1 CYV\ (In(z^) = 0). High stress occurs predominantly over large 

'roughness elements', but only under certain conditions as yet undetermined. 

2 

The Cd-^ = (U./U-^) vs ln(z ) plot is very closely fit by the profile 

lU ^ 10 o 

relationship, 

(using K = 0,40; Businger, et al, 1971, have recently noted that ^ = 0.35 
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is a better value of Von Karman’s constant for neutral stability over a 
land surface); this is merely a check on the consistency of the logarith¬ 
mic profile assumption and on the continuity of data from different time 
periods, and serves to emphasize the significance of the anomalies of 
the other plots (i,e., they appear not to be due to instrumental or analysis 
peculiarities, and the 0,94 correlation is a sufficient criterion for 
assuming a logarithmic distribution). The relationship 

shows that turbulent production increases with increasing fluctuations, but 
with a very strong bimodality, indicative of the presence of more that one 
production-generation regime. 

Fig. 5 shows JH(Uj^q, In(z^) ) and JH(U^q, Cd^^g) for averaging periods 
of 2 and 10 minutes. The basic patterns of anomalies at specific values 
of persist at the longer averaging times, but with less variance. (The 
10-minute averages have been compared by Groscup [1971] to those of Ruggles 
[1970] [on the basis of loglO (z^)]). The basic patterns remain, irrespec¬ 
tive of averaging time, but the decreasing intensity of the anomalies shows 
that, as in the cat’s paw studies, the strong and intermittent events are 
blurred or obscured by long time averaging. Even with these long averages 
and highly logarithmic profiles, it remains clear that the ’mean’ is a 
very poor characterization of physical reality. The time history of events 
and the regimes and instantaneous characteristics of wind and wave fields, 
are highly significant in determining the nature of the boundary layer dy¬ 
namics. 

This point is further emphasized in Fig. 6, which shows JH of with 
U^, In(z^), and Cd^^^ for runs 109 - 159, and JH(Uj^q, U^) for all runs for 
the two subclasses | ^ doA <r »o cm/^cc. The JH of 

runs 109 - 159, characterized by strong and increasing wind and a large wave 
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field, occupy concentrated segments of the overall diagrams. Signifi¬ 
cantly, most of the anomalies are missing, and as previously noted Cd^^ 
for these runs was 0.9 x 10 lower than the overall average; the value 
for these runs is in rough agreement with that* found by most investiga¬ 
tors. This set of runs apparently satisfies the criteria for a single 
or homogeneous wave generation regime (the bimodalities of the other rela¬ 
tionships also disappear or decrease when these runs are plotted separate¬ 
ly) . The division by mean velocity derivative magnitude shows that the 
parameter relationships do vary with time history, and that (as in the 
cat's paw study) higher wind speeds are in general associated with stron¬ 
ger fluctuations. 

Fig. 7 shows examples of the results of the JH of the parameters of 
the logarithmic underwater profile (runs 109 - 159). 30-second mean 
current speed (at lOm) was approximately 10 cm/sec, with a range of -3 
to 27 cm/sec. The diagrams show quite consistent relationships among 
the profile values, with no indication of the anomalies or bimodalities 
found in the atmosphere. values are negative, as the wind drift is 

strongest near the surface. Joint histograms between air and water pro¬ 
file parameters from this set of runs failed to show any overall sta¬ 
tistical relationships; this is not surprising in view of the fact that 
most of the momentum from the wind goes directly into the wave field 
(Stewart, 1967; Dobson, 1971), so that the effect on the surface drift 
current is not immediate or else (see section B.2, Figs. 28 and 29) causes 
the short average profiles to become highly nonlogarithmic. 

In summary, the mean velocity profiles in both air and water are 
predominantly logarithmic during conditions of neutral stability. Ano¬ 
malies in the atmospheric profile parameters demonstrate the importance of 
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time history and the existence of more than one wind-wave generation 
regime, with enhanced activity at certain critical wind speeds under as 
yet unknown circumstances. Anomalies stand out most strongly and the 
air boundary layer becomes much less logarithmic at short averaging 
times, indicating the importance of short, intermittent events in the 
development of unstable shear layers and associated instabilities. The 
profile can not be well represented by a single similarity parameter, and 
there is no simple relationship between wind velocity and the profile 
measures. Attempted use of such relationships for characterizing the 
wind regime and predicting wave growth can lead to significant errors. 
Consideration must be given to the time history of generation and to the 
importance of intermittency, i.e., to the enhanced air-sea coupling and 
momentum transfer at certain wind speeds or in short events. 
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B. Time History 

The 51, 10 minute/hour runs (consecutive except for 1 hour gaps 
between runs 20 and 21 and 42 and 43) selected for detailed analysis were 
taken from 21 to 23 October, 1970. The wind direction had switched from 
N to S 30 hours prior to the first run, as a large high pressure system 
(> 1040mb at the center) passed the site moving NE. Winds increased and 
weather deteriorated throughout the three day period, then winds slacked 
and weather cleared on October 24. Analysis of current records and of the 
tide record at Wood’s Hole, Massachusetts, showed that the surface drift 
continued NE, into Vineyard Sound, after the NE flood tide on the morning 
of October 21. The buoy thus faced into the wind and waves throughout 
this set of runs. This was confirmed by visual observation. 

B.1. 5 and 10 Minute Averages 

Figs. 8 and 9 depict the time changes of the 10 minute 
average horizontal wind speed at 10m (computed) and 2.4m (measured) and 
the wind extrema at 2.4m. The figures also show current speed at Z = 
-1.52m, and stress ( r ryo ) computed from air and water profiles. Short 

lines connect the averages of 2, 5 minute periods from each run and give 
an indication of the variability on shorter averaging times (the large 
offsets in Runs 111 and 142 are due to noise; those in 120 and 157 are 
real; 2.4m velocity extrema from runs 129-133 were unavailable). 

There are strong periodicities in all of the variables. Mean wind 
shows an overall increase from 8-12m/sec; winds increased slowly through 
run 132, then increased rapidly to lOra/sec, stayed near that level for 10 
hours, decreased and remained moderate for a similar period, and then 
again increased rapidly starting around run 152. Superimposed upon this 
long modulation is a shorter but strong periodicity of roughly 2-4 hours. 
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Velocity shear and extrema also show an overall increase with time, with 
superimposed shorter periodicities. The water current has a strong semi¬ 
diurnal tidal component with alternating times between maximum flood of 
5 and 7 hours; the current maxima around run 139 is strongly suppressed. 
Some of the wind fluctuations occur with the tidal frequency and its 2nd 
harmonic; local wind maxima frequently occur during maximum flood. 

The effect of the tidal modulation is vividly displayed in the 
water boundary layer profile as an excess during ebb, when the wind and 
tidal components are opposing and relative wind to water speeds a maximum. 
The short 2-4 hour fluctuations of the wave field are also apparent in the 
water stress. Air stress has a diurnal periodicity and a small 2-4 hour 
component as well as local maxima during the two major wind increases; it 
also temporarily increases in runs 141 - 146 when the winds decrease ra- 
pidly. Mean stresses are thus largest when there is a significant distur¬ 
bance of the equilibrium balance between current, wind and wave fields. 

The stress values from air and water are in rough agreement. Differences 
may be partly accounted for by the momentum carried by propagating swell 
and wind waves, and by the momentum and equilibrium inbalances caused 
by the tides. 

Mollo-Christensen (personal communication) has noted that except for 
the tidal variations, 10 minute mean air and water velocities are compar¬ 
able when normalized by their dynamic viscosities, i.e., . 

This is shown'in Fig. 10 (from Mollo-Christensen and Groscup [1971]). 

This relationship can be accounted for by assuming that the main shear is 
in the viscous sublayer, so that 
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If ^ Sa. > ^•^•9 if sublayer thickness in both media is controlled 

by the same elements (say the wavelets), then since the stresses are in 
rough balance the boundary layers are shear coupled and related by the 
ratio of their dynamic viscosities, k)a ' * The long term ave- 

rage values thus again support the assumption that the major air sea coup- 
ling can be simply described, in spite of the known intermittent discon- 



Figure 10 

term variations of air and water stress. 

Fig. 11, Fig. 12 and Fig. 13 show the variability with time of the 
wave field parameters. Fig. 14 shows the corresponding wave spectra, and 
Fig. 15 the development of selected frequency components (cross sections 
of Fig. 14). 80% and 90% confidence limits at ^ = 10^ are shown by ticks 
inside the ^ axis of Fig. 15. E is mean energy density, E = 
and momentum density M = E/c (c = wave celerity) computed by spectral 
summation. c ) and M ) are respectively the celerity of the 

wave of maximum spectral density and the momentum computed by assuming 
that the energy density is concentrated at that one frequency. The 
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Energy Density Growth Rote, and 
Celerity of wave of Spectral Maximum 
Normalized with Friction Velocity 


Figure 13. 


Runs 109-159 
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information in these figures may be summarized as follows. 

During the two-day period the mean energy density of the wave field 

increased by nearly an order of magnitude, from 2 x 10^ to 17 x 10^ ergs/ 

2 

cm . The period of the wave of maximum spectral density (T) shifted to 
lower frequencies, from 4 seconds at run 109 to ^ 8 seconds at the end. 
These changes can be clearly seen in Fig. 11, Fig. 12 and Fig. 14. The 
changes in T and E are slow and comparatively steady for the first 23 
hours. E then increased rapidly in conjunction with the rapid wind in¬ 
crease that started around run 132. Fig. 14 shows that the spectral peak 
turned sharply toward lower frequencies during this velocity increase; 
after run 132 the separation of the spectrum into two peaks (swell and 
locally generated waves) becomes more distinct and the peaks continue 
to separate as the swell peak shifts to still lower frequencies. The 
equilibrium range above the local peaks has a very rough f ^ slope 
(Phillips [1958]) with occasional local maxima in the range of 0.5 to 1 
hz. The distinct separation between local and swell peaks was wiped 
out during the rapid wind increase of runs 148 and 149. 

All of the figures show that the growth of the wave field, in terms 
of both mean energy density and of individual spectral components, is not 
steady. There are large, significant variations of E, T, normalized par¬ 
ameters, and of )> spectral density at selected frequencies 

(Fig. 15). E (1 Hour) (Fig. 13) becomes as large as + .1% of the hourly 
average E during later runs. The hourly changes, both growth and decay, 
are largest when the v/ave field is most energetic and during major wind 
speed changes. 

Fig. 15 demonstrates that even the higher frequency components (.6310 
and .3550 hz), which are always above the spectral peak and in the ’equili- 
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brium’ or saturation range, have significant (at the 90% confidence level) 
increases and decreases of energy content on these one hour averaging 
periods. The .2515 hz component, which was at the spectral peak at the 
start of the runs, appears to ’overshoot’ in runs 110-112. Its spectral 
density then falls off to a lower, ’equilibrium’ value during the slow 
growth period that extends through run 132. After that run, the 4 second 
waves again show major growth, ’overshoot’ slightly at the same spectral 
density they had earlier reached, and then settle into a second, higher 
’equilibrium’ level. In runs 154 - 156 their energy content increases to 
still higher levels. The wave field can apparently support more than one 
overshoot and saturation range under different wind conditions. 

The .1583 and .1254 hz components both grow until around run 139, 
when the higher of the two frequencies reaches a quasi-steady state level. 

The most rapid growth of the low frequency waves occurs in the rapid wind 
increase starting at run 132. The major wind decrease of runs 142 - 143 
is not accompanied by significant decay; the decrease, and E(1 hour) 

for this period, are no larger than those found in the minor fluctuations 
during growth periods. The fluctuations of the individual frequency 
components is not in phase; some portions of the spectrum grow, others 
decay, in what appears to be a random fashion except during the major 
wind increase at run 132, which is followed around run 135 by growth of 
all but the highest frequency components. The locally generated waves, 

.2515 and .3550 hz, show similar periodicities to those of the local wind 
field. The individual swell and high frequency components are essential- 
independent of these short term variations, while the total wave field - 
E, T and M - does show the effect of local generation. 
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A comparison of the plot of ) with those of E and M 

shows that spectral development is most rapid when the wind is moving fas¬ 
ter than the waves of the spectral peak. Some growth still occurs even 
when the waves move faster than the wind (where Miles theory is clearly 
inapplicable due to the absence of a ’critical layer’). The general 
range of ) is quite low for the entire set of runs. These 

features indicate that swell was propagating into the area, but do not ac¬ 
count for the lack of local decay of low frequency components between runs 
137 and 158 when ) was continuously less than 1. It may be 

hypothesized that energy was being transferred to the low frequency waves 
by nonlinear resonant interactions within the wave field, (e.g., Garrett 
[1970]) or that wind generation occurs at these frequencies during brief 
but intense periods when U was larger than C )• 

The plot of C ( ^ max^^^* shows a strong diurnal periodicity, and 

lags the similar variation in X -by roughly 1^/2. There is no ap- 

a rr 

parent relationship between this parameter and (5 between 

it and energy density growth rates. 

The normalized wave momentum plot indicates that only 60-80% of the 

total wave momentum is carried by waves at the spectral peak; most of the 

remainder is carried by waves at higher frequencies (lower phase speeds). 

Fig. 16 shows the momentum spectrum (5 minute averages, 6 dof) from several 

runs (ordinates have been normalized to the spectral maximum). These 

plots show the overall increase in M and M , the shift of M to lower 

max max 

frequencies, the bi- and trimodality associated with the separation of 
swell and local wind wave peaks, and the occasional very strong concentra¬ 
tion of momentum in a single swell frequency. The momentum spectra have 
roughly the same shape as the power spectra, but tend to emphasize the 
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concentration of the wave momentum in one or more narrow frequency bands. 

The 10% momentum levels — the range of f for M(f) .IM — were 
approximately .15 to .55 hz through run 131, and then decreased slowly to 
roughly .09 - .40 hz. Wave momentum was indistinguishable from noise 
level above 1.40 hz near the start of the runs, and this limit decreased, 
starting at run 132, to approximately 1.00 - 1.20 hz. Most of the mo¬ 
mentum spectra had several discrete peaks, and there was usually a consid¬ 
erable difference in the minor features of the 2,5 minute segments in 
each run. Current spectra at -1.5m (not shown) had maxima at both swell 
and locally generated wave frequencies, and were generally coherent with 
both of these wave field components. 

Coherence between horizontal wind speed at 2.4m and waves at the 
spectral peak measured on the 5 minute subintervals) was 

^ 0.65 for 50%, and ^ 0.50 for 97% of the 102, 5 minute samples. 

Phase ( ©orj ) at that peak was TT for 82% of the coherent runs, and 
^ 2 radians for all of them. As noted by Seeshc3tz (1968), wind 

speeds are maximum over the wave troughs. 

Examples of Yuvj and are displayed in Fig. 17. The shift of 

the peak to lower frequencies with time, and its bifurcation into local 
and swell segments, are clearly evident; there was a secondary maximum 
of at the local wave frequency (roughly .32 hz) for about 80% of 

the samples after run 131. In the majority of these cases (note run 
137-2) there was a decrease in Vbetween the two peaks, with a shift 
in from *n to P/ 2 . and then back to V . Winds were generally in¬ 

coherent with the waves ( Xovj statistically indistinguishable from 0) 
above the local wave peak on this averaging time; one notable exception 
is segment 157-1; this run contained a very strong generation event and 
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is discussed in detail in section B.2. In roughly 30% of the runs, was 

fairly high at the D.C. level; in all of these cases, , With the 

exception of the occasional D.C, peak, was low for f < f J 

( ^ max 

Wind spectra did not consistently display distinct peaks at the wave 

^ 4 

spectral maxima (Fig. 18, 80% and 90% confidence limits at = 10 and 
= 10^ are shown on the ^ axis of run 142). They were fairly white 
( < 1 decade dropoff in the .2 to 2 hz band). Nevertheless there was 
nearly always a spectral increase at or around f ^ cases of 

a double-peaked coherence the wind spectra showed a local minimum between 
the peaks. The examples show (f) and uu (f) the same run seg¬ 

ments as in Fig. 17 (they are rough, 6 dof, and the jaggedness is probab¬ 
ly not significant). 

From these facts and figures the following inferences may be drawn. 
Wave generation and spectral development is strongly non-stationary. There 
are significant, rapid changes in spectral growth and these changes are 
poorly represented by long term means; no single parameter or simple com¬ 
bination of parameters can describe the process. Air and water boundary 
layers (both velocities and stress) are viscously shear coupled. Modu¬ 
lations at the tidal frequency and its harmonics ( or of the period 
ATide/ A t) appear in the wind and wave fields; changes in the tidal flow 
disrupt the equilibrium between wind and wave fields and cause modifica¬ 
tion of momentum advection patterns. 

Even though the high frequency portions of the wave field are near 
saturation and account for insignificant portions of the spectral density 
and momentum, they show significant,rapid (hourly) variations and can 
account for spectral growth when the relationship would normally 

predict major continuous damping. Measured growth rates (long term 
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averages) of these high frequency components are inherently incapable of 
displaying major changes, since on the mean these waves are near satura¬ 
tion, and also because the very high frequencies are continuously being 
generated by the larger waves (Longuet-Higgins,[1963];Crapper [1970]). 

The importance of these high frequency waves lies in their ability to 
couple the wind with larger waves and to assist the growth of low fre¬ 
quency instabilities by nonlinear coupling between scales. It will be 
shown in the next section that these components have extremely rapid 
growth rates in intermittent generation * events’, but again these impor¬ 
tant physical momentum transfer features are obscured by the averaging 
process. 

The influence of waves on the turbulent atmospheric boundary layer 
is strong at both the swell and the locally generated spectral peaks. 
Furthermore, the 2-4 hour periodicity of wind and both the locally gen¬ 
erated and overall wave fields are strongly related. Changes in measures 
of the wind and wave fields (e.g., E (1 hour)) are largest when the mag¬ 
nitude of the measures is largest, similar in nature to the pattern of 
higher frequency components during strong gusts that was noted in the 
cat’s paw study. 

The major conclusion from this section is that study of the time his¬ 
tory of generation and spectral development is essential for an understan¬ 
ding of the physics involved. Wind, wave and current fields show evidence 
of strong interrelationships that are poorly represented by parameteriza¬ 
tion of the mean. The rapid variations in E, E (1 hour), and t all 

indicate that generation is predominantly an intermittent process, char¬ 
acterized by rapid transfer of momentum and followed by weaker couplings 
and nonlinear interchange between spectral components. 
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Figure 21. 
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^•2. 30 Second Averages . 

Several of the 10 minute runs were chosen for de-- 
tailed analysis to help determine the nature and scales of intermittency 
and the resultant momentum transfer. Analysis and data display were simi¬ 
lar to those used for the 10 minute runs. 

Xig. IQ ^nd F1g. _2Q depict the changes with time of 30 s econd 
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Figure 22 

averages of variables in the air and water for runs 110, 127, 157 and 159. 

Fig. 21 shows E (30 sec) and C ($ )/U. for the same runs. The major 

max ^ 

features of the 10 minute runs also stand out in these displays of short 
term change. There is a general increase in U, E, M, T, wave height and 
spectral density and momentum extrema. Stress computed from the air and 
water profiles are in rough agreement, although they show different perio¬ 
dicities. The short term changes are largest and most rapid when the mag¬ 
nitudes are largest; changes in U are paralleled by increases in E, though 
not as clearly as in the longer averages because of the influence of wave 
groups (these had a roughly 40 second periodicity for later runs, and also 
significantly affected the measurements). again normally over 

the troughs, but with much scatter; this is illustrated in Fig. 22, JH 
(U« ., Y^) from 2 short segments, and is more clearly shown in the coherence 
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diagrams (Figs. 29 and 30). 

In these runs, the waves of the spectral peak carry roughly 80% of 
the total momentum density (again implying that a significant portion is 
carried by waves of higher frequency, or slower phase speed), while the 
waves of the momentum peak (when different from the spectral peak) over¬ 
account for the total momentum. U^^/C(6 ) is in the range of 1 - 2, 

10 max 

and is in fact largest during the earliest run, due to the absence of 
swell. This parameter shows no significant peaks at the points of wind 
and E maxima; this may be interpreted to mean that the wind is strongest 
over the largest wave components, indication of selective coupling. 

C ($ xnax^^^* shows no distinct relation to spectral growth or decay; 

E (30 sec) has comparable variations in runs 110 and 127, and the para¬ 
meter is similar in magnitude in runs 110 and 159 when the energy growth 
patterns are radically different. 

For these short averaging times E(30 sec) is 1% to 10% of the mean 
energy density for the corresponding 10 minute periods. The rapid changes 
in spectral shape and content are well shown in Fig. 23, from run 157. 
Cross-sections of this diagram are. shown in Fig. 24 (n.b. 80% and 90% 
confidence limits at 10^ on axis). It is immediately apparent that 

the lower frequencies (referring here both to swell [.1254 hz] and to 
local wind generated waves [.2515 hz]) are relatively constant, with 
variations primarily attributable to propagating wave groups (this was 
checked by examination of analog strip chart displays of the filtered 
signals). The change in the high frequency components, however, is 
very dramatic, with 30-second variations of up to 2 orders of magnitude 
(3 orders of magnitude for the 4.46 hz waves. Fig. 23). The change of 
shape of the ’equilibrium range’ in Fig. 23 also demonstrates this 
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variation; there is a repeated pattern of large excesses in specific high 
frequencies followed by a filling out of the spectrum. The subjective im¬ 
pression of high frequency excesses climbing the slope of the spectrum 
and transferring their energy to lower frequency components is hard to avoid. 
Fig. 24 also shows that the excesses do not alvjays occur simultaneously in 
all frequency bands and that on the 30 second averaging period the high fre¬ 
quency growth does not always correspond to the presence of low frequency 
peaks. Further, intermittent high growth rate events are of short duration, 
typically lasting between 30 and 60 seconds; i.e., similar to the period 
of the major wave groups. 

Changes in the wave field and the wind-wave relationships throu^ 
an 'event* are illustrated for the burst of segment 8, run 157, in Figs. 

25 and 26. Fig. 25 shows the coherence of U and rj from before, during, 

t\ 

and after the event (segments 6-10); Fig. 26 shows the momentum spectra 
from the same period. Before the event is high and =-TT for the 

low frequency components; is insignificant for most higher frequencies. 

During the event Yuy| increases dramatically in the'1.25-3.25 hz range, 
shows a continuous 'slip' through.this range; this is probably accounted 
for by the averaging time being much larger than the high frequency insta¬ 
bilities so that both generation and detritus is represented, and by the 
deceleration of the gust parcel during generation. Immediately after the 
event Voyj has decreased for frequencies >1 hz;*it is still lower in segment 
10, except for the local peak at 1.25 hz where the waves lead the wind (an 
opposite phase relationship from that of segment 9). Before the event the 
wave momentum is strongly concentrated in the local waves (0.32 hz) and 
sv/ell (0.15 hz). During the burst the high frequency components, out to 
over 3 hz, carry a much larger portion of the momentum than at any other 
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Figure 23 
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Figure 24. 
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Figure 25. 
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time. Immediately after the event the very high frequency components have 
disappeared, and the 0.45 hz component has been strongly augmented. This 
peak shifts to a still lower frequency CO.40 hz) in the following segment, 
and eventually nearly all the momentum is again concentrated in the local 
waves and swell. 

Fig. 20 and Fig. 21 show that the event, on the 30 second averaging 
basis, was notable by anomalously high values of wind shear,; wave height, 

current maximum, E, E(30 sec), and M. T increased slightly during the 
event, and significantly 60 seconds later. ^ is not excessively large, 
but the profile for segment 8 had a low correlation with the model 
(= 0.91) so that the computed value of T is unrealistic. The large 
departure of profiles from the logarithmic mean, and the development of 
unstable shear layers, are characteristic of such intermittent events (as 
noted in the cat’s paw study and further shown in Figs. 28 and 29). 

Fig. 27 depicts the detailed (analog strip chart output) structure 
of such a burst type generation event from run 116. Isotachs (from profile 
measurements), wind velocity fluctuations at two heights (high pass filter¬ 
ed at 0.05 hz), wave height, and a measure of the energy content in low and 
high frequency bands of wind and waves are depicted. Scales are relative. 

The major burst occurs after a rapid wind increase and is associated with 
the coupled occurance of low frequency wave groups in air and water. The 
generation of atmospheric turbulence Cvery high frequency wiiid fluctuations) 
occurs predominantly over the wave peaks in air and water, and these peaks 
occur simultaneously. Both low frequency and high frequency wind and wave 
energy are maximum during the burst, and some intermittent activity con¬ 
tinues after the major gust had passed. The low frequency wave components 
are slightly damped after the gust leading edge, when the burst is strong- 
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Figure 27. 
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Figure 29. 
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est. This figure clearly demonstrates a very strong coupling between 
media and between widely separated frequency components. 

Figs. 28 and 29 are similar diagrams from the first and third E peaks 
(see Fig. 20) of run 159 (the ^ scale is comparable in Figs. 27 and 29; 
all others are relative). These figures again show that intermittent, high 
growth rate wave and turbulence generation occurs during the coupled 
occurrence of large, low frequency groups in wind and surface wave fields. 
The event is preceeded by profile distortions, a wind decrease then rapid 
increase. In Fig. 28 the low frequency components are again seen to be 
slightly damped during the intense high frequency instabilities. Large, 
low frequency components of the wind and wave fields that occur indepen¬ 
dently do not lead to this intense coupling and momentum transfer to high 
frequency components; it is the strong nonlinear coupling between the 
media which results in the generation of high frequency instabilities and 
enhanced coupling of scales. 

The strong shear layers in the atmospheric profiles are of very 
short duration. They are almost invariable associated with large wave 
groups and high frequency wave energy spikes, and then quickly revert to 
logarithmic form as the mean wind energy is drained and the atmospheric 
dynamics recover from the local instabilities. The current profiles show 
the effects of the momentum of the large wave groups; they have little 
direct relationship to the instantaneous wind profiles. 

The coherence and phase diagrams for the last event of run 159 are 
shown in Figs. 30 and 31. Most of the large energy peak for this time 
period V7as due to the presence of an exceedingly large wave group, and 
the high frequency generation was less intense than in the similar set from 
run 157. Nevertheless, the sequence of events associated vj^ith the instabil- 
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ity are again apparent. There is the same pattern of high at low 

frequencies and momentum carried in long waves, followed by enhanced high 
frequency coupling and temporary momentum excesses in higher frequencies. 

These figures have demonstrated the short duration, strong coupling, 
rapid growth rate, highly nonlinear and intermittent nature of wave gen¬ 
eration and spectral development. Even the 30 second averages blur the 
intensity of the coupling by overlapping between generation events and 
detritus or propagation. The microstructure is of critical importance, 
both because the majority of momentum transfer occurs at high frequencies 
and because the high frequency instabilities and waves serve as catalysts 
in the coupling and low frequency interactions between scales and media. 

Long sequences of measurements and long averages can predict the background 
conditions under which generation occurs, and can show the effects of very 
low frequency coupling, but are incapable of depicting the essential 
physics of wind-wave generation. 

Summary of Results and Conclusions 

The cat’s paw study showed the.importance of ’bursts* in initial 
generation of high frequency waves and turbulence and indicated the signif¬ 
icance of intermittency and strong nonlinear coupling. The data from 
this study, where a much larger dynamic range of scales is involved, 
support those conclusions and the following statements. 

1. Wave generation is intermittent and involves strong non¬ 
linear coupling among a wide range of frequency components in the wind and 
wave fields. The high frequency instabilities are critical to the momen¬ 
tum transfer process, and much of the air to sea momentum flux goes directly 
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into hxgh frequency waves. The important scales are those of the major 
waves and wave groups, and those of the high frequency wavelets coupled to 
them. These facts are substantiated by the large and rapid fluctuations 
in E, M, and E for both 1 hour and 30 second averaging periods (Figs. 11-13, 
19, 20, and 23), by the simultaneous occurrence of large magnitude, coher¬ 
ent fluctuations of wave height, velocity, and energy content in widely 
separated frequency bands in air and water (Figs. 25, 27-30), and by the 
resultant excesses in spectral density and momentum in high frequency 
waves (Figs. 23, 24, 26, 31). 

2. The high frequency wavelets act as drag elements to couple 
the mean wind to the major wave field; they also aid the growth of lower 
frequency instabilities by nonlinear stress-strain coupling. This was 
demonstrated in t^ie cat’s paw data, and is shown herein by the variations 

in coherence patterns between wind and waves (Figs. 25, 30) and the coupled 
occurrence of bursts in widely separated frequency bands (Figs. 27-29). The 
results of Wu, who makes a similar contention, are not supported by the 
short terra profile averages, primarily because of the complicated dynamics 
of air sea interaction. 

3. The wave field is modulated by surface currents; this 
perturbation in turn causes an imbalance in the equilibrium between wind 
and wave fields. The wind and surface current‘fields are viscously shear 
coupled and their profiles yield comparable mean stress‘values (Figs. 8-10). 

4. There appear to be a variety of wave generation regimes. 
These regimes can be only poorly represented by mean parameters or by 
simple fits to the mean profile conditions. The different regimes are 
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characterized by bimodal distributions of profile parameters (Figs. 2-5) 
that are not present when a ’homogeneous' set of runs is considered separ¬ 
ately (Fig. 6). Profile parameters for a non-homogeneous set of runs have 
larger variances and a much more complicated pattern of values than do 
those of a homogeneous set, and there are no simple formulas which can 
adequately describe the relationships between them. 

5. Study of time history is important for an understanding 
of spectral development. The time history is poorly represented by mean 
parameters or long averages. The approach of the profile to a logarithmic 
form with increasing averaging time obscures the importance of the insta¬ 
bilities associated with short, intermittent perturbations and development 
of shear layers (Figs. 28, 29). The long average spectra (Figs. 14-16) 
can depict the general growth of the wave field but fail to reveal the 
significant flux to high frequencies during intermittent generation (Figs. 
23, 24). The rough quantitative agreement between stresses and viscously 
normalized velocities in air and water boundary layers is strongly modulated 
by tidal forces. 


6. Boundary layer profiles in air and water are essentially 
logarithmic for averaging times > 30 seconds and for neutral stability. 

The profile parameters and Zo can not be expressed in terms of a single 
parameter. The criterion^ 0.94 for a fit of.measured values to the model 
is satisfactory (Fig. 4). There do exist anomalous or ’critical’ wind 
speeds, in the sense that around these speeds the variability in profile 
structure is increased and the resulting profile parameters have anomalous 
means, variances, and amplitude distributions (Figs. 3-5). The mechanics 
responsible for these anomalies have not yet been demonstrated. 
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7. The growth of the major portion of the spectrum is hy¬ 

pothesized to be due both to the wind-wave coupling enhanced by the high 
frequency wavelets and to nonlinear transfer of energy within the wave 
field. The shift from high to low frequency content is best illustrated 
by the momentum spectral changes in Figs. 26 and 31. Spectral growth can 
occur even at ^ (Figs. 11, 12) and this fact is unexplained 

by either propagation or by available theories. 

8. Both the overshoot phenomenon and the equilibrium range 
are meaningful only in terms of long averages. The wave field can support 

a range of equilibrium levels under different wind conditions (Fig 15, par¬ 
ticularly the .2515 hz component), and the wave field is never in a steady 
saturation state (Figs. 15, 23, 24). 

On the basis of these results it is concluded that a large portion 
of wave generation occurs intermittently and is due to strong coupling 
between a large range of scales of motion and between wind and waves. Local 
instabilities and direct nonlinear interactions between waves and turbu¬ 
lence are crucial factoiS in generation, and most of the momentum flux 
from these ’bursts’ goes directly into high frequency waves. 
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IV. CONCLUSIONS AND^ RECOMMENDATIONS 


The questions posed in the introduction may be answered as follows: 

1. The important time scales in the transfer of momentum from 
wind to waves are those of and within the intermittent turbulent bursts. 
These in turn are determined by nonlinear coupling to the period of the 
major waves (localized instabilities or subbursts) and the wave groups 
(the gusts or bursts themselves). The transfer is rapid, produces high 
growth rates of high frequency waves, and quickly drains the available 
energy from the mean wind field. Nonlinear interactions, and transfer of 
energy within the v/ave field by breaking and dissipation, operate on 
similar time scales. Transfer is effectively completed within two or 
three periods of^the waves or wave groups. The time scale of equilibrium 
imbalance and restoration between the mean wind, wave, and current fields 
is on the order of 2 to 4 hours. 

2. The fluctuations involved in wind wave generation are 
intermittent and strongly nonlinearly coupled, involving interactions 
between a wide range of scales. The short instability fluctuations, and 
the resultant energy density growth rates of high frequency waves, are 
roughly three orders of magnitude larger than those observed on the mean. 
The generation of waves and the generation of turbulence are inseparable. 

3. There are several types of coupling between media. These 
include the strong interactions during instability fluctuations; the form 
drag of the wavelets on the mean wind flow and the resultant coupling 
between low frequency, large amplitude wave groups; the interaction of 
wave induced Reynolds stresses at swell and locally generated wave fre- 
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quencies with turbulent Reynolds stresses; an undetermined anomalous 
interplay between wind and waves at discrete critical wind speeds; and 
viscous shear coupling between mean wind and surface drift. Many of these 
interactions appear to depend strongly on the presence and position of 
high frequency fluctuations (the microscales). 

Coupling between scales involves both the weakly nonlinear resonant 
wave interactions, and the strongly nonlinear stress-strain relationships 
involved in the interactions between Reynolds stress and shear. These 
latter operate both within and between media. High frequency instabilities 
aid the growth of lower frequencies; fluctuation magnitudes are largest 
when the range of scales involved in the fluxes are largest. 

4. Generation ’events’ may be distinguished by the short 

time average app::oach to, and deviations from, lognormal distribution of 

velocity fluctuations; by strong coherence between high frequency waves 

and turbulent Reynolds stress or horizontal velocity; by temporary 

excesses of energy and momentum in high frequency (> 1 hz) surface waves; 

and by the enhancement in magnitude and frequency of stress and buoyancy 

fluctuations caused by the strong nonlinear couplings. These features can 

not as yet be parameterized in terms of the mean or of profile parameters, 

but can be used for conditional sampling in a real time data collection/ 

2 

analysis system. The criterion T ^ 100 dynes/cm , when smoothed over 
a time period.^ 0.5 seconds, is a good first approximation for a switch 
to determine the presence of a burst within a cat’s paw. 

5. There are no easily measurable similarity parameters 
which can describe the dynamics of air-sea coupling. The apparent effect 
of fluctuations depends strongly on the averaging scale over which they 
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are evaluated, i.e., the effect is different for different physical 
scales. There are more than one distinct regimes of generation and there 
are large anomalies in the relationships between wind field parameters. 

There is a gross, qualitative similarity in the appearance of weak and 
strong fluctuations, and of large bursts and subbursts, but the intense 
nonlinear interactions and resulting enhancement of the magnitude of 
fluctuations and fluxes preclude any strict similarity representation of 
the physics of momentum transfer. 

The mechanics of the development of a wind-wave field can be heuris- 
tically described in terms of these features: 

Initial generation is caused by the forcing of random turbulent 
pressure fluctuations (the Phillips resonance mechanism); the critical 
time scales involved in this coupling are those of the waves generated. 

In the presence of a strong thermal stratification the generation of waves 
is delayed until the boundary layer is destabilized. The work done on the 
waves creates reaction stresses in the air which generate small scale 
turbulence. The small scale fluctuations do work on the larger scale 
shear in air and water, which aids the development of larger scale 
instabilities and stronger fluctuations. The energy from the mean flow 
is rapidly depleted through transfer to the water surface and dissipation 
to heat. The subbursts occur predominantly within larger Scale anomalies 
or gusts that are initially triggered by atmospheric dynamics. 

After an isolated gust has passed, buoyancy forces rapidly reestablish 
boundary layer stability. The packet of free surface waves generated by the 
gust and bursts appears as roughness elements to the wind. In the presence 
of increasing wind or recurrent gusts the disruption of flow caused by the 
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roughness elements acts to enhance the probability of occurrence of 
instabilities over the packet. The scales of interaction between media 
have increased to include not only the bursts and subbursts but also the 
group or patch of roughness elements and the major waves within it. In 
successive interactions the scales are similarly increased and the 
energy in the initial waves and groups which enhanced the coupling is 
augmented. The scales are yet further enlarged by nonlinear resonant wave 
interactions within and among the wave groups. The very high frequency wave¬ 
lets, which received much of the momentum during the localized high fre¬ 
quency instabilities, either break or otherwise transfer their energy 
to other components, or are ’organized’ by the convergences and divergences 
of orbital velocity in the underlying longer waves. The concentration of 
high frequency components on the crests and leading faces localizes the 
drag elements in a favorable position, since stresses exerted upon them 
are in the same direction as the orbital velocity and celerity of the 
longer waves. Under the influence of the drag elements the mean wind 
decelerates over the crests, accelerates over the troughs. 

As the lower frequency waves increase in size their effect on the 
atmospheric boundary layer increases. Fluctuations in the wind field 
become coherent with these components, and wave induced or reaction 
Reynolds stresses develop which both enhance the pressure feedback from 
the mean flow to themselves (Miles’s mechanism) and interact with atmos¬ 
pheric shear and turbulence. The surface.waves are generated and propagate 
in groups, and enhanced perturbations of the boundary layer on the scale 
of the groups causes similar forcing and feedback on that scale. 

The selectively increased forcing of the wind field by the longer 
vzaves and wave groups, and by the roughness elements associated with 
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the major crests, increases the probability that shear layer instabilities 
will be triggered over them. The gusts are thus selectively coupled to 
these areas. During these instabilities, as in a cat’s paw over a calm 
surface, the majority of momentum transfer occurs during the high frequency 
bursts and subbursts, and the momentum flux goes directly into, the high 
frequency wave components. Energy growth rates are locally three orders of 
magnitude larger than the mean, and the available atmospheric energy is 
rapidly depleted. The high frequency waves again break or otherwise 
transfer their energy to the longer waves and swell. The process occurs 
intermittently and can cause spectral growth even when the largest 
waves are moving faster than the mean wind. 

The generation and spectral development process is self-enhancing 
and involves strong coupling between scales of motion in both wind and 
wave fields. During periods of relatively steady winds a quasi-equilibrium 
balance is maintained, but this is easily broken by changes in surface 
current or by wind velocity shifts. A large portion of the momentum 
flux to the surface occurs during the intermittent, large magnitude, 
strongly coupled bursts; and most of this flux goes directly into high 
frequency waves. The microscales play a dominant role both by being the 
active recipients in the momentum transfer and by assisting the growth of 
lower frequency components by stress-strain coupling, by acting as roughness 
or drag elements, and by weakly resonant interactions. 

Many questions concerning the physics of wind-wave generation 
remain unanswered. This study has demonstrated the importance of inter- 
raittency, strong nonlinear coupling, and microscales in momentum transfer, 
but has given no indication of the percent of the total transfer attribu¬ 
table to these features. Further, the effects but not the physics of 
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energy transfer within the wave field have been described. 

Meaningful future experimentation must provide a better description 
of the horizontal and vertical structure of Reynolds stress in air and 
water; quantitative measurement of wave growth during intermittent 
generation and its relation to the development of the wave field; a sta¬ 
tistical description of the structure and spacing of intermittencies and 
their relationship to waves and wave groups; simultaneous measutements of 
high and low frequency wave energy and their spatial and time relation¬ 
ships; a detailed examination of the air-sea coupling at the anomalous 
critical wind speeds; and controlled determination of the effects of 
high frequency degrading slicks. 

Open ocean measurements of velocity and turbulence in the troughs of 
major waves are still unavailable; the role of breaking and the effect 
of the resulting vertical vortices on spectral transfer is poorly under¬ 
stood; and the effects of wind-surface current coupling and the underwater 
boundary layer on wave generation have still been inadequately described. 
Future studies must consider the importance of time history and inter- 
mittency during generation, and most progress will be made through 


detailed, near naviface measurements. 
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APPEiraiX I. INSTRUI'ENTS AND CALIBRATION 

A. Sites and Platforms . 

Buoy Experiment . Instruments vjere mounted on the M.I.T. Air-Sea 
Interaction Buoy (Mollo-Christensen and Dorman,(1971)). The site was at 
the south-western end of Vineyard Sound (4l 20''40" N, 70° 55* 45" W). 

The buoy is shown in Fig. AI.l, The site is 8.9 km south of Cuttyhunk 
Harbour, 7.8 km west of Gay Head Light, Martha's Vineyard, and approxi¬ 
mately 50 km ENE of Block Island. The bottom in the vicinity is very 
fine sand and mud, depth 3I.I ™ S't MLLV/, and mean tidal range less than 
1 m. The bottom is nearly flat to the south and southvjest, and shoals 
slowly to the north and northeast. 

The experiment was conducted from July to November, 1970. Winds in 
the area are dominated by a land-sea breeze regime during the summer, and 
are predominantly southwesterly in the fall. Tidal currents are primarily 
semidiurnal -with a clock-wise ellipse. They are strongly affected by the 
winds. Recorded current maxima -were around 1.5 ra/sec -with a typical slack 
duration of 20 minutes and nearly equal duration of flood (NE) and ebb 
(S^W). 

A shore support and data collection facility was housed in the Coast 
Guard Station (deactivated) on Cuttyhunk. The facility had line-of-sight 
communication -with, and visual observation of the buoy. Ship support was 
provided by the R/v R. R. Shrock (M.I.T.'s oceanographic research vessel) 
and by project owned and local smallcraft. The Shrock was equipped with an 
electronic bathythermograph used to monitor thermal structure during ex¬ 
periments. A small weather station -was maintained on Cuttyhunk. 

The buoy -was slack moored (scope of 2) and consequently faced into the 
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tide (one notable exception being runs 109-159 when continuous high wind 
reversed the surface flow and kept the buoy pointing southwest). Torque 
balance to counter the drag of the underwater instrument tower was accom¬ 
plished by angular preadjustment of the mooring bridle. The bridle was 
placed near the center of lateral resistance so the buoy and attached 
instruments had < 1° mean tilt during all current conditions. 

The buoy and its small motions had little effect on the data. This may 
be shown by the folloiving considerations. Static stability was high due to 
the 4,27 m metacentric height and small surface peretrating cross-section. 

A vertical load of 100 lbs, 10' from the centerline at platform height, 
caused only 0,23° tilt. The large virtual mass and damper plate resulted 
in a critically damped natural heave period of 21,3 seconds, with an am¬ 
plitude response for 6 second waves of A == (Buoy Amplitude)/(Sea Surface 
Amplitude) = 0,074, Yaw was determined by field experience to be of sig¬ 
nificance only during the hour or so at and around slack. Fluctuating 
tilt was measured with a pendulum pot and showed a maximum response at 
around 10 seconds with maximum tilts of 1,5°. Instruments were positioned 
30 ° fon-;ard of the normal to the mean current to minimize interference 
calculated on the basis of potential flovj. Crude estimates on this basis 
indicated that velocity distortion (AU/u) was 0,028 for the minimeters, 
0,009 for the anemometers. The estimates are of course worst near the navi- 
face, in regions of high shear. Wave gauges were placed on both outriggers 
and their outputs compared to eliminate buoy interference during adverse 
wave directions. 

The buoy was equipped with a thermoelectric generator to allow unat¬ 
tended operation for up to 2 weeks under mean condition run schedules 
(12 minutes eash hour), or four hours of continuous operation vjith over- 
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night recharge. Primary power source was a set of eight, 12 volt 120 amp- 
hour nicad batteries. The buoy has been in use for two years and has 
operated successfully in ;'rinds over 20 m/sec. 

Catamaran Experiment , This experiment was conducted in March, 1971j 
from the dock of the M.I.T. Yacht Club boathouse. Bulk atmospheric condi¬ 
tions were generally stable (see Table 1, Section II). Data were useful 
only when winds vjere in the southern quadrants, such that the fetch was 
over 300 

A diagram of the catamaran is shovm in Fig. AI.2. It was designed to 
support up to 1.5 kg of instrumentation in a nearly constant position rela¬ 
tive to the instantaneous surface components of f<1.2 hz. Heave decoupling 
and minimization of mean instrumental tilt were accomplished by the use of 
tvjo vertical supports, each sliding through tvjo delrin bearings which were 
separated by 0.8 m and affixed to a carefully leveled support structure. 
Pitch and roll bearings vjere located respectively at the attachment point 
of the fore and aft crossbars to the pontoons, and at the vertical support- 
crossbar attachment. This arrangement gave the catamaran a 60 cm vertical 
range and allowed pitch and roll angles of up to 20° with less than 1.5° 
tilt of the sensors. Construction materials were polystyrene foam and .020"^ 
6061 aluminum sheeting and tubing. Total weight was 3.2 kg, draft 5 cm, 
and design natural frequency 1.5 hz. 

Instruments V7ere attached to the forvjard support tube from 15 cm belovr 
to 30 cm above the surface. They were positioned within 5 cm of the bow 
line to minimize flow obstruction. Additional sensors v;ere mounted on the 
support structure. Cables vrere floated aft of the cat to lessen vjeight and 
permit directional control for a full 360 °. Accelerometry tests were con¬ 
ducted v;ith full instrumental complement during most runs. They showed 
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that the cat closely folloxved waves of f1.3 hz, and had negligible re¬ 
sponse to those over 1.85 hz. Waves excited by catamaran motion caused 
extraneous spectral peaks at 50 snd 70 hz but had little effect on frequen¬ 
cies of interest. 

B. Instruments . Most of the instruments were designed and constructed 
at M.I.T. and have seen several years of use and development. In addition 
to precision and accuracy, a major design criterion has been systematiza¬ 
tion of the entire array of transducers to insure that data fomats will 
be compatible and easy to process. Hardware and electronics vdll be des¬ 
cribed in Project reports now in preparation, and also have been discussed 
at length in other reports from the Project (Seesholtz (1968); Ruggles (I 969 ) 
Von Zweck ( 1969 )). The essential characteristics and calibration results 
are as follows. 

Velocity Sensors . Profile velocities were measured -vath a set of 
5 matched cup anemometers in the air and 3 ducted propeller 'minimeters' 
in the water. Both instrument types use chopper v;heels to interrupt a light 
beam that is sensed by photoresistors. Static and dynamic anemometer cali¬ 
brations were conducted in the M.I.T. low turbulence wind tunnel. The 
results indicated that the instruments had a distance constant of 101 cm 
and a static linear calibration of U(cm/sec)=2.325x(pulses/min) + 62 in the 
speed range 0.8 to 30 m/sec. A repeat calibration at the end of the summer 
shovjed deviations from initial values to be less than 2 $, Directional tests 
shovjed that the units had little sensitivity to vertical components. The 
anemometers are dynamically .nonlinear (Wyngaard, et al (1971)) due to a 
faster response to viind increases than to decreases. Overestimates in gusty 
vrinds were roughly 6'j. Miniraeters were calibrated in a tow^ tank at Woods 
Hole Oceanographic Institution, at speeds of 8 to IO 5 cm/sec. The static 
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response was linear, U(cn/sec) = 0.9lllx{pulses/sec) + 2.5. Distance 
constant was 8 cn. Directional response (Fig. AI. 3 ) was best fit by a 
cos© + K sinO curve with K differing between instruments. 

Velocity measurements from the catamaran were made with DISA type 55A81 
industrially coated hot film probes operated as constant temperature ane¬ 
mometers with overheat ratios of 1.55/1 ^.ir and I. 08 /I in water. Static 
calibrations v/ere simultaneously conducted v;ith all 5 probes in the M.I.T. 
low turbulence vrind tunnel. The response vjas expressed in terms of a best 
fit (in the least mean squares sense) quartic polynomial. Directional re¬ 
sponse (Fig. AI. 3 ) was similarly fit by the linear relation 
U(cra/sec) = Q(1.0 - 0.0021(0 -25))> where Q = actual velocity and & - 

attack angle, in the range 25°< © 80° (correlation = 0.92). Dynamic 

response (Rasmussen and Bellhouse (I 968 )) was satisfactory for the frequen¬ 
cies encountered during cat's paws and in fact compensated for the slight 
velocity overestimate at high attack angles inherent in the directional fit. 

Vfave Sensors . Have gauges were capacitance type; they were design¬ 
ed and built at iI.I.T. and have provided satisfactory results in the field 
and laboratory for 4 years, at wave heights up to 5 s^wd frequencies up 
to 100 hz. The output is a variable frequency that is created by beating 
the frequency from an oscillator controlled by the sensor capacitance 
against a standard. The \ 7 ave gauges used on the buoy had a 4.6 m long, .95 
cm diameter stainless steel wire core steering cable sensor, and a 30 .I hz 
standard oscillator. Full range bandvjidth vras 400 hz with resultant re¬ 
sponse of 3.2 hz/cm. The capillary wave gauge for the catamaran had a 15 
cm long, 0.076 cm diameter teflon coated sensor and a 200 hz standard oscil¬ 
lator. Full range variation vjas 1.4 to 2.2 khz, or amplitude response of 
3.88 hz/ram. Sensors i^ere frequently cleansed with alcohol to reduce film 


effects. 
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Tenperature Sensors . Tvro types of thermistors, both v;ith 15 kil 
resistance at 25°C vrere used to pemit interchangability in the bridge 
circuit and provide roughly equal time constants in air and water. Water 
temperature was measured with matched Fenwall K496 oceanographic isocurve 
submini thermistors. Each unit had tvjo, O.O 76 cm diameter bead-in-glass 
tips. Air to still water time constant was 70 ms. The Fenwalls were epoxy 
mounted in I .3 cm diameter stainless steel tubing fitted with retractable, 
ventilated shields. The shields were kept forvjard during the mean condition 
mins (removed during catamaran runs), and decreased the time constant in 
moving water to approximately.2 seconds. Air thermistors on the catamaran 
were VECO model 42A70, .0508 cm diameter bead-in-glass probes mounted in 
.63 cm diameter tapered wooden dowling. Their air to still water time con¬ 
stant vras 35 ws. Using the result of Aagaard (1964) for similar units, 
amplitude response at nyquist frequency (125 hz) for the VECOs was 0.44. 

The differential bridge circuit was designed at M.I.T. It was used in 
the absolute mode vri.th a 6° - 26° C window for mean condition runs; accur¬ 
acy and precision ^'jith both thermistor types was * 0.1°C. Measurements 

from the catamaran used the relative mode (therraistors balanced against 
precision resistors selected on the basis of mean temperature) with a 5^C 
window, and had a precision of ^ 0.02°C. A bomb calorimeter was employed 
as field and calibration temperature standard. 

Additional instruments included a pendulum pot for measurement of buoy 
tilt, a tri-axial strain gauge accelerometer for measurement of catamaran 
motion, and an electronic bathythermograph connected to an x-y plotter 
on the Shrock. 

C. Mounting. 


Buoy . Cup anemometers were attached to aluminum channel which in 
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turn vfas bolted to the meteorological tower. Cups were slightly more than 
1 m from the tower. Heights (for runs 6? - 159) were 2.4, 3*7> 5*0* 7«5j 
and 11.4 m. Ilinimeters were attached to carts vjhich ran on one set of 
tracks on the undervrater tower. They could be depth-positioned by hand from 
the buoy outrigger. Sensors vrere over $0 cm from the tovjer which in turn 
was 2.5 m from the main spar. Mean depths (runs 109 - 159) were -1.52, -2.74, 
and -4.7 m. Thermistors were mounted on the ams of the minimeters and an¬ 
emometers. The primary wave gauge vras located directly above the minimeters; 
the secondary was on the opposing outrigger. Hot film and thermistor elec¬ 
tronics were housed in a vrooden box located between the braces of the 
tower support, and the main electronics package was in a compartment in 
the main spar, directly under the personnel platform. 

Catamaran. The accelerometer and wave gauge electronics were set 
in the fonrard cross-brace. Orthogonal hot films and thermistors v;ere at¬ 
tached to the fon-^ard vertical support •vn.th delrin mounts. The wave gauge 
sensor vras stretched betvreen Wo delrin braces on the forward support. 
Cross-stream separation betvreen the wave gauge and hot film tips was 0.95 
cm; between the vrave gauge (air) or hot films (vrater) and thermistors 1.3 cm 
(see Table 1, Section II). Dovmstream sepatation vras 1 mm. Orthogonal 
positioning of hot films vras insured by pre-fixed mounting and was checked 
by measurements on projected photographs. Separation betvreen hot film tips 
vras 3 ran* 
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APPEIIDIX II. DATA COLLECTION AND ANALYSIS 

Both hard- and softvjare packages >0.11 be described in detail in forth¬ 
coming reports of the M.I.T. Buoy Project. This section will describe the 
functional nature of the system and outline the major subroutines and alo- 
grithms used in the analysis. The major components of the system are 
depicted in Figs. AII.l and All.2. 

A. Collection and Format . 

Buoy System .. Information from the transducers passes through con¬ 
ditioning amplifiers designed to produce 0-5 volts variable D.C. output. 
This output drives voltage controlled oscillators set on standard IRIG 
channels 3-15 and 17. The channels were matched against anticipated data 
banckri-dth to insure at least a 5^1 bandvjidth ratio. Accelerometer pulses 
(and minimeter pylses during some runs) were sampled at 100 hz and combined 
on a single IRIG channel. A standard 25 khz signal was also transmitted 
and later used for synchronization and pulse separation during breakout. 
Signals from up to 22 separate transducers could be accomodated simultan¬ 
eously. Signals from all IRIG channels were nixed and then used to fre¬ 
quency modulate the transmission signal (243.8 mhz). The major components 
of the buoy data handling pack.age were modified from an Atlas Ilissle Telem¬ 
etry System (USAF excess equipment). 

Data received at the base station on Cuttyhunk was recorded in analog, 
direct mode on a single channel of magnetic tape on a Hewlett Packard 
Model 39603 (KP 396 OB) 4 track recorder. A second.channel was reserved 
for voice narration, weather log, and time signal. Both sides of 2400' 
tapes were used for mean condition runs; one tape held four, 12 minute 
hourly runs. 





Autocorrelation 


160 


05 



Demodulation 










































































































































161 . 


The only analysis attempted at the site was demodulation•and monitoring 
of selected channels to insure signal quality and continuity. The buoy was 
turned on and off remotely by a modified HeathiCit radio control set which 
activated a microsiritch in the buoy electronics compartment. 

Catamaran System . Signals from catamaran transducers were recorded 
on a Princeton Instruments Model 6100 (PI 6 IOO) 8 track recorder. Hot film, 
thermistor, and accelerometer data were passed thi’ough conditioning ampli¬ 
fiers to optimize voltage bias and gain, and then recorded in FIl mode. 
Capillary wave gauge signals, in frequency format, were recorded in direct 
mode. One channel was reserved for voice narration; this vjas subsequently 
used to help select segments for digitization. The incoming data was con¬ 
tinuously monitored to insure signal quality and continuity. 2400* tapes 
allowed over 2 hours of continuous recording per run. 

B. Analysis . 

General . A combination of digital and analog analysis was employed. 
Analog processing with strip chart and x-y recorder display vjas used for 
initial evaluation. Energy content in selected frequency bands was estimated 
by filtering and squaring the signals prior to plotting. Analog spectral 
analysis of potentially interesting events was performed with the assistance 
of a PI correlator and Fourier analyzer and the x-y plotter. The initial 
processing familiarized the user with the data and permitted early estimates 
of fluctuation magnitude for later use in digital processing. 

Three digital computers were used for secondary and final analysis. 
Digitization, preprocessing, and simple analyses such as pulse counts and 
calibration checks were performed on the Project’s Honej’^jell H -316 compu¬ 
ter. This unit has an 8 k memory and is equipped vrith a I 6 channel Raytheon 
’miniverter* analog to digital (a/d) converter Model MADC 12A, The miniver- 
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ter outputs 12 bit a/d words for a full range of i 10 volts;-least bit pre¬ 
cision is 5 wv. Programming was in Honeytrell's DAP-16 Assembler level 
language. Output from the Honeywell v/as in two forms: the digitized data 
in bit-string, a/d word format on 7 track magnetic tape, and pulse count 
and calibration information on either the ASP.-33 teletypei'Triter or paper 
tape. The latter information was transferred to punched cards. 

A PDP-7 computer (at Linclon Laboratory) equipped I'ri.'th railtiple tape 
drives and cathode ray tube (CRT) display was used to examine in detail 
the tapes produced by the 11-316. The bit-string data was first converted to 
geophysically meaningful units; statistics were calculated; values from or¬ 
thogonally mounted hot films vrere converted to U and ¥; and the resulting 
information was then rerecorded in a format standard to the facility. The 
resultant data were displayed on the CRT, limited filtering and spectral 
analysis performed, and hard copy graphical output obtained for further 
study. Programming was done in Fortran IV by Mr. B. Laird. 

Final analysis vras done on the IBM 36 O at the M.I.T. Inforraation Pro¬ 
cessing Facility (IPC). Programs were written in PL/i to reduce l/O prob¬ 
lems, and extensive use was made of a Stromberg-Carlson 4020 (SC-4020) plot¬ 
ter for display of time series, probability distributions, and spectra. 

The data handling system is capable of statistical, time history, and 
frequency domain analysis of up to I 6 channels of input data simultaneously. 
Analog preprocessing capabilities include rerecording at a choice of speeds 
and formats, filtering, integration, squaring, orthogonal decomposition, 
amplification, and pulse separation, as vrell as totally analog computation 
and display of time history, correlations, and auto-and cross-spectra. The 
H -316 allows digitization at up to 250 hz. 

Buoy Data . The mixed signals from the HP3960B were separated in a 
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set of 8 demodulators (from the Atlas Ilissle System). Sampled anemometer 
signals vjere separated using ring counters and the 25 khz clock signal for 
timing. Sections of the data vrere filtered, squared, amplified and display¬ 
ed on a 4-track Sanborn recorder to provide an initial overview. Tito sep¬ 
arate digitizing passes vrere then made. The first yielded velocity profile 
infonnation from a pulse counting program on the H-3l6. Tvrenty, 30 secon d 
counts from 10 minutes of each mean condition run were made, transferred to 
punched cards, and run through a two part program on the IBM '^ 60 » The first 
digital program computed profile parameters on the basis of the logar¬ 
ithmic profile model. The second program used punched card output from the 
first to compute single and joint histograms and statistical information, 
displayed on SC-4020 plots. Velocity values from selected runs (vjith aver¬ 
aging times as short as 5 seconds) vjere hand plotted in profile, isotach, 
and vertical shear format for comparison with the strip chart analog data. 

In the second digitizing pass, the runs selected for time history anal¬ 
ysis v?ere rerecorded on the PI 6 IOO in demodulated form. The data vjas fil¬ 
tered to avoid aliasing (velocity values at 2.5 hz and v;ave height at 5 hz 
using matched Roclajell Model 1010? operational filters), amjDlified to the 
full range of the miniverter, and digitized at 10 hz. Calibration during 
this preprocessing was maintained by passing sinusoidal and D.C. signals 
through the entire system. The calibration data vTOre displayed on the ASR-33 
and transferred to cards for use in the main analysis. Additional calibra¬ 
tion checks during digital processing were made by comparison of pulse 
count and digitized velocity values. The time signal initially recorded on 
the voice channel of the HP39^>OB vjas used to synchronize the tvjo passes. 
Final analysis used the softv;are package described belovr. 

Catamaran Data . All of the cat’s paw data vjas pre-vievred using 
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analog techniques. Selected portions were then digitized at 250 hz; the 
data had been low-pass filtered at 100 hz using the internal PI 6 IOO filters. 

The voice narration was used to aid in distinguishing the presence of gusts 
until their inprint on the analog signals became readily recognizable. 

Nine tapes containing some 2.5 x lo"^ data points from each of 8 channels 
were obtained. These tapes contain infonnation from over 75 well identified 
cat’s paws as ^^ell as a large amount of background data. The tapes ^^ere t 

examined on the L.L. PDP-7 and segments selected for final analysis. 

IHl 360 Software. The analysis package was designed to permit 
rapid, standardized analysis of time series information sampled at equal 
intervals. Flexibility and user control were major design criteria. The 
package can handle up to I 6 channels of information simultaneously, and 
allows for the inclusion of two additional channels of computed values such 
as Reynolds stress and buoyancy fluctuations. In addition, the initial 
information can be replaced during analysis by other computed values. The 
package may be of use to other investigators, and is consequently briefly 
described herein. 

Program flovj is diagrammed in Fig. All.2. It is highly modular and 
designed to handle a wide variety of time series information, with a mini¬ 
mal number of routines specifically restricted to air-sea interaction 
studies. Flexibility is maintained by user option control from punched card 
parameters. There is provision for a user supplied routine (OTHER) which 
can be reached from both time and frequency domains. Time domain analysis 
is performed first; in the frequency domain segment the original data 
(stored in a three-dimensional matrix) is replaced by its fourier transform 
coefficients. The major features and alogrithms are as follows; details iri.ll 
be given in a forthcoming Project report. 
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COMTRL handles tape positioning, dimensions the major data matrix, 
selects plot/no plot mode, and inputs data from the tape. This is the only 
routine that is specifically linked to the H -316 7 track magnetic tape 
format. 

DOIT is called by CONTRL upon completion of input. This subroutine 
does little computing, but determines vjhich other routines vrill be called 
and the sequence of their use. It does contain the banning and stress for¬ 
mation routines. 

REA.LVAL performs, with operator option, the following functions: 
conversion of a/d words into voltage; conversion from voltage to real values 
with up to a quartic fit (required for linearization of the hot films); de¬ 
termination of trend, extrema, and the 1st 4 central moments; computation 
and plotting of a histogram with user-supplied number of intervals and 
limits set either to computed extrema or user-supplied; adjustment of values 
to a calibration standard (used for comparison of pulse count and digital 
velocities) and/or removal of mean (supplied or computed) and/or trend. 

Ilean and trend may also be removed in DOIT immediately prior to the fre¬ 
quency domain segmaht. 

FILTER operates by fast fourier transforming the data, modifying or 
deleting the selected frequency components, and then retransforming. The 
user supplies the modification factors so that effectively any desired win¬ 
dow may be employed. A standard cosine taper of adjustable length is inclu¬ 
ded. This routine can be called from up to 3 different locations in DOIT 
so that the data may be used in both filtered and unfiltered forra. 

ORTHOG computes horizontal (U) and vertical (VJ) velocities from ortho¬ 
gonally mounted instruments. The decomposition is done in taro steps; first 
an estimate of U and VI is obtained from the known directional response and 
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velocity readings of the tuo instruments, and then an alignment correction 
is made to compensate for mean tilt or oscillatory fluctuations of the 
probe holder. The background and details of this procedure were discussed 
at length by Ruggles (1969). There is a choice of either a cos0 + K sin© 
or a linear directional fit to the data. The alignment correction assumes 
that V/ = 0 over a number of data points selected by the user, and that the 
instruments are truly orthogonally mounted. A known mean horizontal vel¬ 
ocity, if available, can be used for additional correction. The choice of 
the averaging interval for V/ = 0 is critical, and is made by a considera¬ 
tion of the dynamics of the system and by post-facto examination of plots 
of the data to insure continuity of values at interval end points. For the 
catamaran data, 60 points xjere used for air hot films, 250 for water 
hot films. 

PROFIL computes parameters of the assumed logarithmic profile model; 
effectively, it compares U(z) to (U*/k) IuCz/zq) by a best fit in the least 
mean squares sense, taking the instmiment heights as known (In previous 
vjork from this project, e.g. Ruggles (1969)) velocity values were taken as 
knovm, z as the variable. In the initial stages of this analysis both methods 
were used, and for profiles sho'.’ing correlations with the model better than 
0.9^ the resulting differences in U* and Zq were always less than 5/^*). 
Reynolds stress values are available, the production term of the horizontal 
component of turbulent kinetic energy equation (-u'w* dU/dz) can be calcu¬ 
lated. The separate program to handle the pulse count data uses similar 
alogrithms. 

JOIiITP computes and plots joint histograms, and also provides one-way 
histograms and means-on-interval of the input variables. 

TIIIEPLT plots variables against time, with user controlled or computed 
scaling and number of strips per SC-4020 frame. It can be accessed both 
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before and after filtering. This subroutine also plots histograins. 

F0URT2 is a USASI Basic Fortran version of the Cooley-Tukoy fast fourier 
transfora, vrritten by II. Brenner. It can handle real or complex data in 
multidimensional, arbitrarily sized arrays. End values are dropped and 
tapered in DOIT to provide a ■vrindovi prior to transforming; in this analysis 
approximately lO^j of the data vjas so treated to limit Gibbs effect. The 
returned coefficients are normalized and banned in DOIT, then a density 
factor is applied to result in the power spectral density estimate. The 
banned estimate has a distribution id-th 6 dof; restriction of averaging 
to 3 values vras necessary in the catamaran study to preserve detailed 
resolution in view of the high sampling rate and relatively few points 
available for analysis (core and averaging limitations). 

XSPECTR utilizes the noraalized fovirier transform coefficients to 
determine cross-spectral densities (Co and Quad spectra). These eu'e then 
used to compute coherence magnitude ( ) and phase (© ). The averaging 

interval used to compute coherence is selectable by the user; in this study 
9 points were used for samples under 500 points, 15 for the 5-roinute mean 
condition run segments, and 25 for the 10-second catamaran time periods. 

The results were then banned for the longer samples. 

SPECOUT prints and plots spectral data. Autospectra were plotted on a 
log(magnitude) linear (frequency) grid (log-log is also available). Magni¬ 
tude scaling is based on spectral extrema (up to 10 decades can be accomo¬ 
dated), but the user controls the frequency range. Coherence and phase 
were plotted from the cross-spectral analyses; the user again controls the 
frequency range and number of plots. 

OTHER is the user supplied routine. It can be accessed from two loca¬ 
tions in the time domain and at the end of the frequency domain computations 
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in DOIT. Tivo major versions were used in this study, one each for the buoy 
and catamaran experiments. The buoy version computed mean energy density, 
both directly from wave height data and by spectral accumulation, and wave 
momentum at individual frequency components and by considering the energy 
to be concentrated in waves at the frequency of energy and momentum spec¬ 
tral maxima. The catamaran version had four segments; one computed buoyancy 
fluctuations and determined joint histograms between buoyancy and other 
variables; a second formed the natural logarithm of squared and normalized 
velocity derivatives for the lognormality study; a third created the switch 
function that was used to evaluate intermittency; and the fourth contained 
a band pass filter for rapid examination of filtered signals. 

Computational Kethods . Let xi, i = 1, 2, ... N represent the 
sampled data values where HAt = T is the total sample length and l/At =A f 
the sampling frequency. Then for any a/d word yi, voltage 
Vj_ = -10.0 + yj^ X 204.8 from the knovm characteristics of the miniverter, 
and Xi = % Aj v^ in geophysical units where the -ij's are the conver¬ 
sion values from calibration. Then 


mean = ^ T 


i ^ 

j^^ moment = 

i 

H 

2 (<; - x’V SO that 

; s i ' 

variance = 


— X _ 

- x* = crj- 

^ *> 

skewness = 


- 3x5 x B 5^ , and 

kurtosis = 

IT? 

— y + X — s • 


Trend (T) is found from a best fit in the least mean squares sense to the 
data. Define Vj sTl i ^ , and minimize 

->'r) - (if; -T i- -/sj^ ^ where both Xj_ and i are considered 

to be random variates. Then the expectation of 5 is 
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E(s^W E(;<y;_-o-Tr-n - (/9-x 

= o;'^+‘l!V> -^(|3-v -r-TTTf- xTcs^k; 

vrhere is the covariance of and i. The,positive definite part must 

= 0 for minimum S , thus /S = x - X i . 

C'xi 

V/ith the correlation coefficient ^ 

adding and subtracting /ocr^, E(S*') = {X - pis^f *(^-P'^)g^ 


so that for minimum ^ , T = ~ and {i-P)g],'^ 

M ah 

— ;v - — 

N fii *■ " '■ * 

Thus, 


T 


H ^ 


but 

so that 


i. = - > /Z. I - 


N(.K>-v>')(N + i) 
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and if the mean and trend are removed, 


5; - CHr' -)) 


For the calculation of profile parameters, the assumed model is 

u,... 

where U* is friction velocity, K. is Von Karman's constant (= 0.40) and 
Zq is roughness length. Using the least squares method developed above, 
substituting Uj_ for x^, S = U#/kl for T* , ln(z) for i and 
(-Uy./K. ) In(zo) for /S , 


e 


^ *^.4 (s') 


S = Jw(0 
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Further, stress "c = DUj): V7here D= density of the medium, 

U(ic.m\ r \J,„ = ^3 X 2.LOIS’S5 - 

for velocity in cm/sec and heights in cm, ahd drag coefficient 

Figure All,3 sho^^s the assumed mounting and notation for orthogonal 

decomposition. Assuming that the hot 
films have a linear directional response 
in some appropriate range of , let the 
velocity recorded be Vi, so that 

V7here Qi is the actual velocity and P 
and f are determined from calibrations. 

For the DISA hot films, =25°, P = .0021 
tfs and the relationship holds {/° = 0,92) 

in the attack angle range 25 °^ ^ 80°, 

But 



- 

45“ . 0 

and 

=45“ - 0 

SO that 


Q. (l_T>(-2-C. - 6;)) 


= (l- 7 (zo-»- 0.)) ^ 


Q. = 

y.; ^6-. 

> 0.' 




Z{|- zot) 

T(w ''VO 



Then horizontal (U) and vertical (V/) velocity estimates are obtained from 

V/i = si*4 e, ^ Q;Cose; . 

For minimeters, the directional response is best -given by 
V; - (co-s -V '7'si.i <i>C) 

in the range 0° < ^ 80°, xfith the value of P varying betaijeen instru¬ 

ments, Substituting for ^ in terras of 6' and then in terns of W and U, 

r {I y. . -K ~ 

‘ z V , + -p J ^ “'t-Vi-?;' 
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To correct for tilt, assume 17^^ and found as above are actually at 
an angle c<; to the correct cLxis. Then essentially, letting and'V repre¬ 
sent the corrected values, 

*Vj i = v/c Co ■So/' - ^ PiCe.-5cW; 

Require VI = 0; then if a known U is available, the correction to lowest 
order is found from 


I 


<\J > 
u 


while if U is not knowiq tam 


Co^o^; = 

C VJ> 
“ 4 0> 


±}y> 

u 


is used. The averag¬ 


ing times are selected by the user to minimize mean or fluctuating tilt. 
Total rms error in computed Reynolds stress using this lox^est order tech¬ 
nique can be as large as 20'^. For very small values of , say < 10°, 
the linear hot film directional response fit overestimates 'I± by roughly 
1.5'^j* Assuming io'j nis turbulence levels for norraal atmospheric conditions, 
the probable rms deviation of around 45° is roughly 10°; thus this type 
of fit overestimates only extremely high turbulence values. This is somevjhat 
compensated for by the frequency response of the- units, but still tends 
to emphasize anomalously high stress values which is useful for detection. 

In spectral computations, F0URT2 returns the real and imaginary coef¬ 
ficients of the discretized fourier transform, 

•>--72 • • 

iThere t = j ht and s = -feAf, (" i /S’ is returned. 

The coefficients are normalized by the number of data points, e.g. 
aC'b'i = z«;(- 5 ') , and the power or variance spectral estimate is formed: 


^ u') -- iF«r - dwN . 

X >y 
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Hanning is applied: 

M. u 


a. ^ * -L (4'> * a 

A -V. ■' 1. ^1<« A 


and the power spectral density estimate is fonued by multiplying by the 
density factor, d6-t -T 


i.e. 




where is the Nyquist or folding frequency = l/2^t. 

Computation of the cross-spectrum (follovjing Wunsch, (1970).) is based 
on the hypothesis that the variables are related by the convolution 


relationship, 


2 


A- 


i-iz 


for all realizations of x and y, taking x and y as random variates so that 
^ is a knovm transient. Then 

where ID and G are the fourier transforms of % and y. Then 

For the hypothesized convolution relationship to have meaning, 6(s) must 
be a smoothly varying function of s,'and coherence is a measure of this 
smoothness, i.e., ^j^ 


C<.v,w ' 


Z ■=' 


<0 f'S*"') 


For peaked spectra this is best computed in a form which includes a weighing 
by the power at the various frequencies, i.e.. 




1F. ts) ] 


v?here denotes the complex conjugate. But 
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Fx U) r I , and jf. (s)] iG^Cb^l -- ( 


Coherence nas^nitude is thus defined as 






Using p and i to denote real and imaginary parts of the fourier trans¬ 
form, the alogrithm is 

z */ 

r M L + ('F’l (?»«■)<5^(5+*^)^] *" 




Similarly, phase 


6 l ‘>) -- ■: 7 t^n-' (-2— 


M is chosen by the user and end values are assumed symmetrical for the 
averaging. The alogrithm is designed so that positive means 

y leads x. 

Mean wave energy density is computed from the relationship 

t - 

vihere yj is wave height, /“ density, and <3 gravity. But = ^77 
where ^ is the autocorrelation function (not computed); in turn, 

0 ^ 6 ^ OC> ^ 

-o. ~U> -<• 

Both direct and spectral summation techniques vrere used. The results vrere 
alTfays vjithin 1 . 5 ^ of each other for spectra xjith 6 dof or more. Departures 
of 20^5 vrith mean departure roughly 5 ^j occured with 2 dof (raw or unhanned 
short spectral estimates). 
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Ilonentun ;:as computed from the relationship 

v/here c is linear wave celerity, 

C'(s') = ^TaAl^Uh') 5,' Ta.U 

Jk = vjave number and h = 31 Thus, 

-- /’S §«,7^^Vcu) 

Mean momentum density iras computed both by summation over the spectral 
components and by assumption that E was primarily attributable to the wave 
of maximum variance or momentum, i.e,, 

M- - -i- 
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